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METHODS AND APPARATUS FOR 
PREDICTING, CONFIRMING, AND DISPLAYING 
FUNCTIONAL INFORMATION DERIVED FROM GENOMIC SEQUENCE 



CROSS REFERENCE TO RELATED APPLICATIONS 

5 The present application: is a continuation- 

in-part of U.S. patent application serial nos. 
09/632,366, filed August 3, 2000 and 09/608,408, filed 
June 30, 2000; claims the benefit under 35 U.S.C. 
§ 119(e) of U.S. provisional patent application serial 

10 nos. 60/236,359, filed September 27, 2000, 60/234,687, 
filed September 21, 2000, 60/207,456, filed May 26, 
2000, and 60/180,312, filed February 4, 2000; and 
further claims the benefit under 35 U.S.C. § 119(a) of 
GB patent application no. 0024263.6, filed October 4, 

15 2000, the disclosures of which are incorporated herein 
by reference in their entireties. 
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FIELD OF THE INVENTION 

The present invention is in the fields of 
bioinformatics and molecular biology, and relates 
particularly to analytical methods and apparatus for 
5 predicting, confirming, and displaying functional 
information derived from genomic sequence. The 
invention particularly relates to methods and apparatus 
for identifying portions of genomic sequence data that 
encode genes, to the design, manufacture and use of 
10 genome-derived single-exon nucleic acid microarrays for 
assaying expression thereof, and to methods and 
apparatus for display of genomic sequence annotated 
with expression information. 

BACKGROUND OF THE INVENTION 

15 For almost two decades following the 

invention of general techniques for nucleic acid 
sequencing, Sanger et al., Proc. Natl. Acad. Sci. USA 
70 (4) : 1209-13 (1973); Gilbert et al . , Proc. Natl. Acad. 
Sci. USA 70 (12) : 3581-4 (1973), these techniques were 

20 used principally as tools to further the understanding 
of proteins — known or suspected — about which a basic 
foundation of biologic knowledge had already been 
built. In many cases, the cloning effort that preceded 
sequence identification had been both informed and 

25 directed by that antecedent biological understanding. 

For example, the cloning of the T cell 
receptor for antigen was predicated upon its known or 
suspected cell type-specific expression, by its 
suspected membrane association, and by the predicted 

30 assembly of its gene via T cell-specific somatic 
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recombination. Hedrick et al . , Nature 308 (5955) : 149-53 
(1984) . Subsequent sequencing efforts at once 
confirmed and extended understanding of this family of 
proteins. Hedrick et al . , Nature 308 ( 5955) : 153-8 
5 (1984) . 

More recently, however, the development of 
high throughput sequencing methods and devices, in 
concert with large public and private undertakings to 
sequence the human and other genomes, has altered this 

10 investigational paradigm: today, sequence information 
often precedes understanding of the basic biology of 
the encoded protein product. 

One of the approaches to large-scale 
sequencing is predicated upon the proposition that 

15 expressed sequences - that is, those accessible through 
isolation of mRNA - are of greatest initial interest. 
This "expressed sequence tag" ("EST") approach has 
already yielded vast amounts of sequence data. Adams 
et al., Science 252:1651 (1991); Williamson, Drug 

20 Discov. Today 4:115 (1999); Strausberg et al . , Nature 
Genet. 15:415 (1997); Adams et al . , Nature 
377 (suppl . ) : 3 (1995); Marra et al . , Nature Genet. 
21:191 (1999). For nucleic acids sequenced by this 
approach, often the only biologic information that is 

25 known a priori with any certainty is the likelihood of 
biologic expression itself. By virtue of the species 
and tissue from which the mRNA had originally been 
obtained, most such sequences are also annotated with 
the identity of the species and at least one tissue in 

30 which expression appears likely. 

More recently, the pace of genomic sequencing 
has accelerated dramatically. When genomic DNA serves 
as the initial substrate for sequencing efforts, 
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expression cannot be presumed; often the only a priori 
biologic information about the sequence includes the 
species and chromosome (and perhaps chromosomal map 
location) of origin. 
5 With the ever-accelerating pace of sequence 

accumulation by directed, EST, and genomic sequencing 
approaches - and in particular, with the accumulation 
of sequence information from multiple genera, from 
multiple species within genera, and from multiple 
10 individuals within a species - there is an increasing 
need for methods that rapidly and effectively permit 
the functions of nucleic sequences to be elucidated. 
And as such functional information accumulates, there 
is a further need for methods of storing such 
15 functional information in meaningful and useful 

relationship to the sequence itself; that is, there is 
an increasing need for means and apparatus for 
annotating raw sequence data with known or predicted 
functional information . 
2 0 Although the increase in the pace of genomic 

sequencing is due in large part to technological 
changes in sequencing strategies and instrumentation, 
Service, Science 2 8 0 : 995 (1998); Pennisi, Science 283: 
1822-1823 (1999), there is an important functional 
25 motivation as well. 

While it was understood that the EST approach 
would rarely be able to yield sequence information 
about the noncoding portions of the genome, it now also 
appears the EST approach is capable of capturing only a 
30 fraction of a genome's actual expression complexity. 

For example, when the C. elegans genome was 
fully sequenced, gene prediction algorithms identified 
over 19,000 potential genes, of which only 7,000 had 
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been found by EST sequencing. C. elegans Sequencing 
Consortium, Science 282:2012 (1998). Analogously, the 
recently completed sequence of chromosome 2 of 
Arabidopsis predicts over 4000 genes, Lin et al. r 
5 Nature, 402:761 (1999), of which only about 6% had 

previously been identified via EST sequencing efforts. 
Although the human genome has the greatest depth of EST 
coverage, it is still woefully short of surrendering 
all of its genes. One recent estimate suggests that 
10 the human genome contains more than 146,000 genes, 

which would at this point leave greater than half of 
the genes undiscovered. It is now predicted that many 
genes, perhaps 20 to 50%, will only be found by genomic 
sequencing . 

15 There is, therefore, a need for methods that 

permit the functional regions of genomic sequence - and 
most importantly, but not exclusively, regions that 
function to encode genes — to be identified. 

Much of the coding sequence of the human 

2 0 genome is not homologous to known genes, making 
detection of open reading frames ("ORFs") and 
predictions of gene function difficult. Computational 
methods exist for predicting coding regions in 
eukaryotic genomes. Gene prediction programs such as 

25 GRAIL and GRAIL II, Uberbacher et al . , Proc. Natl. 

Acad. Sci. USA 88 (24) : 11261-5 (1991); Xu et al . , Genet. 
Eng. 16:241-53 (1994); Uberbacher et al . , Methods 
Enzymol. 266:259-81 (1996); GENEFINDER, Solovyev et 
al., Nucl. Acids. Res. 22:5156-63 (1994); Solovyev et 

30 al. r Ismb 5:294-302 (1997); and GENE SCAN, Burge et al . , 
J. Mol. Biol. 268:78-94 (1997), predict many putative 
genes without known homology or function. Such 
programs are known, however, to give high false 
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positive rates. Burset et al . , Genomics 34:353-367 
(1996) . Using a consensus obtained by a plurality of 
such programs is known to increase the reliability of 
calling exons from genomic sequence. Ansari-Lari et 
5 al., Genome Res. 8(l):29-40 (1998). 

Identification of functional genes from 
genomic data remains, however, an imperfect art. For 
example, in reporting the full sequence of human 
chromosome 21, the Chromosome 21 Mapping and Sequencing 
10 Consortium reports that prior bioinf ormatic estimates 
of human gene number may need to be revised 
substantially downwards. Nature 405:311-199 (2000); 
Reeves, Nature 405:283-284 (2000). 

Thus, there is a need for methods and 
15 apparatus that permit the functions of the regions 

identified bioinf ormatically - and specifically, that 
permit the expression of regions predicted to encode 
protein - readily to be confirmed experimentally. 

Recently, the development of nucleic acid 
2 0 microarrays has made possible the automated and highly 
parallel measurement of gene expression. .Reviewed in 
Schena (ed.), DNA Microarrays : A Pra ctical Approach 
(Practical Approach Series ) , Oxford University Press 
(1999) (ISBN: 0199637768); Nature Genet. 
25 21 (1) (suppl) : 1 - 60 (1999); Schena (ed.), Microarrav 
Biochip: Tools and Technology , Eaton Publishing 
Company/BioTechniques Books Division (2000) (ISBN: 
1881299376), the disclosures of which are incorporated 
herein by reference in their entireties. 
30 It is common for microarrays to be derived 

from cDNA/EST libraries, either from those previously 
described in the literature, such as those from the 
I.M.A.G.E. consortium, Lennon et al . , "The I.M.A.G.E. 
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Consortium: an Integrated Molecular Analysis of Genomes 
and Their Expression, Genomics 33(l):151-2 (1996), or 
from the construction of "problem specific" libraries 
targeted at a particular biological question, R.S. 
5 Thomas et al . , Cancer Res. (in press). Such 

microarrays by definition can measure expression only 
of those genes found in EST libraries, and thus have 
not been useful as probes for genes discovered solely 
by genomic sequencing. 
10 The utility of using whole genome nucleic 

acid microarrays to answer certain biologic questions 
has been demonstrated for the yeast Saccharomyces 
cerevisiae. De Risi et al - , Science 278:680 (1997). 
The vast majority of yeast nuclear genes, approximately 
15 95% however, are single exon genes, i.e., lack introns, 
Lopez et al., RNA 5:1135-1137 (1999); Goffeau et al., 
Science 274:563-67 (1996), permitting coding regions 
more readily to be identified. Whole genome nucleic 
acid microarrays have not generally been used to probe 
20 gene expression from more complex eukaryotic genomes, 
and in particular from those averaging more than one 
intron per gene. 



SUMMARY OF THE INVENTION 

25 The present invention solves these and other 

problems in the art by providing methods and apparatus 
for predicting, confirming, and displaying functional 
information derived from genomic sequence. 

In one aspect, the invention provides a 

30 process for predicting functional regions from genomic 
sequence, confirming and characterizing the functional 
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activity of such regions experimentally, and then 
associating and displaying the information so obtained 
in meaningful and useful relationship to the original 
sequence data. 

5 In a related aspect, the present invention 

provides apparatus for verifying the expression of 
putative genes identified within genomic sequence. In 
particular, the invention provides novel genome-derived 
single exon nucleic acid microarrays useful for 

10 verifying the expression of putative genes identified 
within genomic sequence. 

In another aspect, the present invention 
provides compositions and kits for the ready production 
of nucleic acids identical in sequence to, or 

15 substantially identical in sequence to, probes on the 
genome-derived single exon microarrays of the present 
invention . 

In further aspect, the present invention 
provides a genome-derived single-exon microarray 

20 packaged together with such an ordered set of 

amplifiable probes corresponding to the probes, or one 
or more subsets of probes, thereon. In alternative 
embodiments, the ordered set of amplifiable probes is 
packaged separately from the genome-derived single exon 

25 microarray. 

In another aspect, the invention provides 
means for displaying annotated sequence, and in 
particular, for displaying sequence annotated according 
to the methods and apparatus of the present invention. 

30 Further, such display can be used as a preferred 

graphical user interface for electronic search, query, 
and analysis of such annotated sequence. 

In another aspect, the invention provides 
genome-derived single exon nucleic acid probes useful 
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for gene expression analysis, and particularly for gene 
expression analysis by microarray. The invention 
particularly provides genome-derived single-exon 
probes known to be expressed in one or more tissues. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other objects and advantages of 
the present invention will be apparent upon 
consideration of the following detailed description 
taken in conjunction with the accompanying drawings, in 

10 which like characters refer to like parts throughout, 
and in which: 

FIG. 1 illustrates a process for predicting 
functional regions from genomic sequence, confirming 
the functional activity of such regions experimentally, 

15 and associating and displaying the data so obtained in 
meaningful and useful relationship to the original 
sequence data, according to the present invention; 

FIG. 2 further elaborates that portion of the 
process schematized in FIG. 1 for predicting functional 

20 regions from genomic sequence, according to the present 
invention; 

FIG. 3 illustrates a visual display according 
to the present invention, herein denominated a 
"Mondrian", in which a single genomic sequence is 

25 annotated with predicted and experimentally confirmed 
functional information; 

FIG. 4 presents a Mondrian of a hypothetical 
annotated genomic sequence, further identifying typical 
color conventions when the Mondrian is used to annotate 

30 genomic sequence with exon-specif ic expression data, as 
in FIGS. 9 and 10; 
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FIG. 5 is a chart that summarizes data from 
experimental Example 1, showing the size distributions 
of predicted exon length (dashed line) and actual PCR 
products (amplicons) (solid line) as obtained from 
5 human genomic sequence according to the methods of the 
present invention; 

FIG. 6 is a histogram that summarizes data 
from experimental Examples 1 and 2, showing the number 
of tissues in which predicted exons could be shown to 
10 be expressed using simultaneous two color hybridization 
to a genome-derived single exon microarray of the 
present invention. The graph shows the number of 
sequence-verified products that were either not 
expressed in any of the ten tested tissues/cell types 
15 ("0"), expressed in one or more but not all tested 

tissues ("1" - "9"), or expressed in all tissues tested 
("10") ; 

FIG. 7 is a pictorial representation of data 
from experimental Examples 1 and 2, showing the 

20 expression (ratio relative to control) of probes having 
verified sequences that were expressed with signal 
intensity greater than 3 in at least one tissue, with: 
FIG. 7A showing both the expression as measured by 
microarray hybridization in each of the 10 measured 

25 tissues and the expression as measured 

"bioinformatically" by query of EST, NR and SwissProt 
databases; with FIG. 7B showing the legend for display 
of physical expression (ratio) in FIG. 7A; and with 
FIG. 7C showing the legend for scoring EST hits as 

30 depicted in FIG. 7A; 

FIG. 8 is a chart of data from experimental 
Examples 1 and 2, showing a comparison of normalized 
CY3 signal intensity for arrayed sequences that were 
identical to sequences in existing EST, NR and 
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SwissProt databases (known) or that were dissimilar 
(unknown) , where the dashed line denotes the signal 
intensity for all sequence-verified products with a 
BLAST Expect ("E") value of greater than le-30 
5 (1 x 10" 30 ) ("unknown") and the solid line denotes 

sequence-verified spots with a BLAST expect ("E") value 
of less than le-30 (1 x 10" 30 ) ("known"); 

FIG. 9 presents a Mondrian of BAC AC008172 
(bases 25,000 to 130,000), containing the carbamyl 
10 phosphate synthetase gene (AF154830 . 1) ; and 

FIG. 10 is a Mondrian of BAC A049839. 

DETAILED DESCRIPTION OF THE INVENTION 

Definitions 

As used herein, the term "microarray" and 

15 equivalent phrase "nucleic acid microarray" refer to a 
substrate-bound collection of plural nucleic acids, 
hybridization to each of the plurality of bound nucleic 
acids being separately detectable. The substrate can 
be solid or porous, planar or non-planar, unitary or 

20 distributed. 

As so defined, the term "microarray" and 
phrase "nucleic acid microarray" include all the 
devices so called in Schena (ed*), DNA Microarravs: A 
Practical Approach (Practical Approach Series ) , Oxford 

25 University Press (1999) (ISBN: 0199637768); Nature 
Genet. 21 ( 1 ) (suppl ) : 1 - 60 (1999); and Schena (ed.), 
Microarray Biochip: Tools and Technology , Eaton 
Publishing Company/BioTechniques Books Division (2000) 
(ISBN: 1881299376), the disclosures of which are 

30 incorporated herein by reference in their entireties. 
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As so defined, the term "microarray" and 
phrase "nucleic acid microarray" also include 
substrate-bound collections of plural nucleic acids in 
which the nucleic acids are distributably disposed on a 
5 plurality of beads, rather than on a unitary planar 
substrate, as is described/ inter alia, in Brenner et 
al. f Proc. Natl. Acad. Sci. USA 97 (4) :166501670 (2000), 
the disclosure of which is incorporated herein by 
reference in its entirety; in such case, the term 

10 "microarray" and phrase "nucleic acid microarray" refer 
to the plurality of beads in aggregate. 

As used herein with respect to a nucleic acid 
microarray, the term "probe" refers to the nucleic acid 
that is, or is intended to be, bound to the substrate. 

15 As used herein with respect to solution phase 

hybridization, the term "probe" refers to the nucleic 
acid of known sequence that is, or is intended to be, 
detectably labeled. In either such context, the term 
"target" refers to nucleic acid intended to be bound to 

20 probe by Watson-Crick complementarity. 

As used herein, the expression "probe 
comprising SEQ ID NO", and variants thereof, intends a 
nucleic acid probe, at least a portion of which probe 
has either (i) the sequence directly as given in the 

25 referenced SEQ ID NO, or (ii) a sequence complementary 
to the sequence as given in the referenced SEQ ID NO, 
the choice as between sequence directly as given and 
complement thereof dictated by the requirement that the 
probe be complementary to the desired target. 

30 As used herein, the phrase "expression of a 

probe" and its linguistic variants means that the probe 
hybridizes detectably at high stringency to nucleic 
acids that derive from mRNA. 
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As used herein, the term "exon" refers to a 
nucleic acid sequence bioinf ormatically predicted to 
encode a portion of a natural protein. 

As used herein, the phrase "open reading 
5 frame" and the equivalent acronym "ORF" refer to that 
portion of an exon that can be translated in its 
entirety into a sequence of contiguous amino acids. As 
so defined, an ORF is wholly contained within its 
respective exon and has length, measured in 

10 nucleotides, exactly divisible by 3. As so defined, an 
ORF need not encode the entirety of a natural protein. 

As used herein, the phrase "alternative 
splicing" and its linguistic equivalents includes all 
types of RNA processing that lead to expression of 

15 plural protein isoforms from a single gene; 

accordingly, the phrase "splice variant (s) " and its 
linguistic equivalents embraces mRNAs transcribed from 
a given gene that, however processed, collectively 
encode plural protein isoforms. 

20 For example, and by way of illustration only, 

splice variants can include exon insertions, exon 
extensions, exon truncations, exon deletions, 
alternatives in the 5 T untranslated region ("5 ! UT") 
and alternatives in the 3 1 untranslated region 

25 ("3 ! UT"). Such 3 T alternatives include, for example, 
differences in the site of RNA transcript cleavage and 
site of poly (A) addition. See, e.g., Gautheret et al . , 
Genome Res. 8:524-530 (1998). 

As used herein, the phrase "specific binding 

30 pair" intends a pair of molecules that bind to one 

another with high specificity. Binding pairs typically 
have affinity or avidity of at least 10 7 , preferably at 
least 10 8 , more preferably at least 10 9 liters/mole. 
Nonlimiting examples of specific binding pairs are: 
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antibody and antigen; biotin and avidin; and biotin and 
streptavidin. 

As used herein with respect to the visual 
display of annotated genomic sequence, the term 
5 "rectangle" means any geometric shape that has at least 
a first and a second border, wherein each of the first 
and second borders is capable of mapping uniquely to a 
point of another visual object of the display. 



Methods and Apparatus for Identifying, Confirming, 
10 and Displaying Functional Regions of Genomic 

Sequence 

FIG. 1 is a flow chart illustrating in broad 
outline a first aspect of the present invention, a 
process for predicting functional regions from genomic 

15 sequence, confirming and characterizing the functional 
activity of such regions experimentally, and then 
associating and displaying the information so obtained 
in meaningful and useful relationship to the original 
genomic sequence data. 

20 The initial input into process 10 of the 

present invention is drawn from one or more 
databases 100 containing genomic sequence data. 
Because genomic sequence is usually obtained from 
subgenomic fragments, the sequence data typically will 

25 be stored in a series of records corresponding to these 
subgenomic sequenced fragments. Some fragments will 
have been catenated to form larger contiguous sequences 
("contigs"); others will not. A finite percentage of 
sequence data in the database will typically be 

30 erroneous, consisting inter alia of vector sequence, 

sequence created from aberrant cloning events, sequence 
of artificial polylinkers, and sequence that was 
erroneously read. 
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Each sequence record in database 100 will 
minimally contain as annotation a unique sequence 
identifier (accession number) , and will typically be 
annotated further to identify the date of accession, 
5 species of origin, and depositor. Because database 100 
can contain nongenomic sequence, each sequence will 
typically be annotated further to permit query for 
genomic sequence. Chromosomal origin, optionally with 
map location, can also be present. Data can be, and 

10 over time increasingly will be, further annotated with 
additional information, in part through use of the 
present invention, as described below. Annotation can 
be present within the data records, in information 
external to database 100 and linked to the records 

15 thereto, or through a combination of the two. 

Databases useful as genomic sequence 
database 100 in the present invention include GenBank, 
and particularly include several divisions thereof, 
including the htgs (draft) , NT (nucleotide, command 

20 line), and NR (nonredundant ) divisions. GenBank is 
produced by the National Institutes of Health and is 
maintained by the National Center for Biotechnology 
Information (NCBI) . Databases of genomic sequence from 
species other than human, such as mouse, rat, 

25 Arabidopsis thaliana r C. elegans, C. brigsii, 

Drosophila melanogaster , zebra fish, and other higher 
eukaryotic organisms will also prove useful as genomic 
sequence database 100. 

Genomic sequence obtained by query of genomic 

30 sequence database 100 is then input into one or more 
processes 200 for identification of regions therein 
that are predicted to have a biological function as 
specified by the user. Such functions include, but are 
not limited to, encoding protein, regulating 
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transcription, regulating message transport after 
transcription, regulating message splicing after 
transcription, regulating message degradation after 
transcription, contributing to or controlling 
5 chromosomal somatic recombination, contributing to 
chromosomal stability or movement, contributing to 
allelic exclusion or X chromosome inactivation, and the 
like . 

The particular genomic sequence to be input 

10 into process 200 will depend upon the function for 

which relevant sequence is to be identified as well as 
upon the approach chosen for such identification. 
Process step 200 can be iterated to identify different 
functions within a given genomic region. In such case, 

15 the input often will be different for the several 
iterations . 

Sequences predicted to have the requisite 
function by process 200 are then input into 
process 300, where a subset of the input sequences 

20 suitable for experimental confirmation is identified. 
Experimental confirmation can involve physical and/or 
bioinf ormatic assay. Where the subsequent experimental 
assay is bioinf ormatic, rather than physical, there are 
fewer constraints on the sequences that can be tested, 

25 and in this latter case therefore process 300 can 
output the entirety of the input sequence. 

The subset of sequences output from process 
300 is then used in process 400 for experimental 
verification and characterization of the function 

30 predicted in process 200, which experimental 

verification can, and often will, include both physical 
and bioinf ormatic assay. 

Process 500 annotates the sequence data with 
the functional information obtained in the physical 



and/or bioinf ormatic assays of process 400. Such 
annotation can be done using any technique that 
usefully relates the functional information to the 
sequence, as, for example, by incorporating the 
functional data into the sequence data record itself, 
by linking records in a hierarchical or relational 
database, by linking to external databases, by a 
combination thereof, or by other means well known 
within the database arts. The data can even be 
submitted for incorporation into databases maintained 
by others, such as GenBank, which is maintained by 
NCB1 . 

As further noted in FIG. 1, additional 
annotation can be input into process 500 from external 
sources 600. 

The annotated data is then optionally 
displayed in process 800, either before, concomitantly 
with, or after optional storage 700 on nontransient 
media, such as magnetic disk, optical disc, 
magnetooptical disk, flash memory, or the like. 

FIG* 1 shows that the experimental data 
output from process 400 can be used in each preceding 
step of process 10: e.g., facilitating identification 
of functional sequences in process 200, facilitating 
identification of an experimentally suitable subset 
thereof in process 300, and facilitating creation of 
physical and/or informational substrates for, and 
performance of subsequent assay, of functional 
sequences in process 400. 

Information from each step can be passed 
directly to the succeeding process, or stored in 
permanent or interim form prior to passage to the 
succeeding process. Often, data will be stored after 
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each, or at least a plurality, of such process steps. 
Any or all process steps can be automated. 

FIG. 2 further elaborates the prediction of 
functional sequence within genomic sequence according 

5 to process 200. 

Genomic sequence database 100 is first 
queried 20 for genomic sequence. 

The sequence required to be returned by 
query 20 will depend, in the first instance, upon the 
10 function to be identified. 

For example, genomic sequences that function 
to encode protein can be identified inter alia using 
gene prediction approaches, comparative sequence 
analysis approaches, or combinations of the two. In 
15 gene prediction analysis, sequence from one genome is 

input into process 200 where at least one, preferably a 
plurality, of algorithmic methods are applied to 
identify putative coding regions. In comparative 
sequence analysis, by contrast, corresponding, e.g., 
20 syntenic, sequence from a plurality of sources, 

typically a plurality of species, is input into process 
200, where at least one, possibly a plurality, of 
algorithmic methods are applied to compare the 
sequences and identify regions of least variability. 
25 The exact content of query 20 will also 

depend upon the database queried. For example, if the 
database contains both genomic and nongenomic sequence, 
perhaps derived from multiple species, and the function 
to be predicted is protein coding in human genomic DNA, 
30 the query will accordingly require that the sequence 
returned be genomic and derived from humans. 

Query 20 can also incorporate criteria that 
compel return of sequence that meets operative 
requirements of the subsequent analytical method. 
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Alternatively, or in addition, such operative criteria 
can be enforced in subsequent preprocess step 24. 

For example, if the function sought to be 
identified is protein coding, query 20 can incorporate 
5 criteria that return from genomic sequence database 100 
only those sequences present within contigs 
sufficiently long as to have obviated substantial 
fragmentation of any given exon among a plurality of 
separate sequence fragments. 

10 Such criteria can, for example, consist of a 

required minimal individual genomic sequence fragment 
length, such as 10 kb, more typically 20 kb, 30 kb, 
40kb, and preferably 50 kb or more, as well as an 
optional further or alternative requirement that 

15 sequence from any given clone, such as a bacterial 

artificial chromosome ("BAC")/ be presented in no more 
than a finite maximal number of fragments, such as no 
more than 20 separate pieces, more typically no more 
than 15 fragments, even more typically no more than 

20 about 10 - 12 fragments. 

Our results have shown that genomic sequence 
from bacterial artificial chromosomes (BACs) is 
sufficient for gene prediction analysis according to 
the present invention if the sequence is at least 50 kb 

25 in length, and if additionally the sequence from any 

given BAC is presented in fewer than 15, and preferably 
fewer than 10, fragments. Accordingly, query 20 can 
incorporate a requirement that data accessioned from 
BAC sequencing be in fewer than 15, preferably fewer 

30 than 10, fragments. 

An additional criterion that can be 
incorporated into the query can be the date, or range 
of dates, of sequence accession. Although the process 
has been described above as if genomic sequence 
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database 100 were static, it is of course understood 
that the genomic sequence databases need not be static, 
and indeed are typically updated on a frequent, even 
hourly, basis. Thus, as further described in 
5 experimental Examples 1 and 2, infra , it is possible to 
query the database for newly added sequence, either 
newly added after an absolute date or newly added 
relative to a prior analysis performed using the 
methods and apparatus of the present invention. In 

10 this way, the process herein described can incorporate 
a dynamic, temporal component. 

One utility of such temporal limitation is to 
identify, from newly accessioned genomic sequence, the 
presence of novel genes, particularly those not 

15 previously identified by EST sequencing (or other 

sequencing efforts that are similarly based upon gene 
expression) . As further described in Example 1, such 
an approach has shown that newly accessioned human 
genomic sequence, when analyzed for sequences that 

20 function to encode protein, readily identifies genes 
that are novel over those in existing EST and other 
expression databases. In fact, as shown below, fully 
2/3 of genes identified in newly accessioned human 
genomic sequence have not hitherto been identified. 

25 This makes the methods of the present invention 
extremely powerful gene discovery tools. 

And as would be appreciated, such gene 
discovery can be performed using genomic sequence from 
species other than human. Particularly useful species 

30 are those used as model systems during drug 

development, such as rodent, particularly mouse. 

If query 20 incorporates multiple criteria, 
such as above-described, the multiple criteria can be 
performed as a series of separate queries or as a 



single query, depending in part upon the query 
language, the complexity of the query, and other 
considerations well known in the database arts. 

If query 20 returns no genomic sequence 
meeting the query criteria, the negative result can be 
reported by process 22, and process 200 (and indeed, 
entire process 10) ended 23, as shown. Alternatively, 
or in addition to report and termination of the initial 
inquiry, a new query 20 can be generated that takes 
into account the initial negative result. 

When query 20 returns sequence meeting the 
query criteria, the returned sequence is then passed to 
optional preprocessing 24, suitable and specific for 
the desired analytical approach and the particular 
analytical methods thereof to be used in process 25. 

Preprocessing 24 can include processes 
suitable for many approaches and methods thereof, as 
well as processes specifically suited for the intended 
subsequent analysis . 

Preprocessing 24 suitable for most approaches 
and methods will include elimination of sequence 
irrelevant to, or that would interfere with, the 
subsequent analysis. Such sequence includes repetitive 
sequence, such as Alu repeats and LINE elements, vector 
sequence, artificial sequence, such as artificial 
polylinkers, and the like. Such removal can readily be 
performed by identification and subsequent masking of 
the undesired sequence. 

Identification can be effected by comparing 
the genomic sequence returned by query 2 0 with public 
or private databases containing known repetitive 
sequence, vector sequence, artificial sequence, and 
other artif actual sequence. Such comparison can 
readily be done using programs well known in the art, 
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such as CROSS_MATCH or RE PEATMASKER , the latter 
available on-line at 

http : //ftp . genome . Washington . edu/RM/RepeatMasker . html , 
or by proprietary sequence comparison programs the 
5 engineering of which is well within the skill in the 
art . 

Alternatively, or in addition, undesirable, 
including artifactual, sequence can be identified 
algorithmically without comparison to external 

10 databases and thereafter removed. For example, 

synthetic polylinker sequence can be identified by an 
algorithm that identifies a significantly higher than 
average density of known restriction sites- As another 
example, vector sequence can be identified by 

15 algorithms that identify nucleotide or codon usage at 
variance with that of the bulk of the genomic sequence. 

Once identified, undesired sequence can be 
removed. Removal can usefully be done by masking the 
undesired sequence as, for example, by converting the 

20 specific nucleotide references to one that is 
unrecognized by the subsequent bioinf ormatic 
algorithms, such as "X". Alternatively, but at present 
less preferred, the undesired sequence can be excised 
from the returned genomic sequence, leaving gaps. 

25 Preprocessing 24 can further include 

selection from among duplicative sequences of that one 
sequence of highest quality. Higher quality can be 
measured as a lower percentage of, fewest number of, or 
least densely clustered occurrence of ambiguous 

30 nucleotides, defined as those nucleotides that are 
identified in the genomic sequence using symbols 
indicating ambiguity. Higher quality can also or 
alternatively be valued by presence in the longest 
contig. 
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Preprocessing 24 can, and often will, also 
include formatting of the data as specifically 
appropriate for passage to the analytical algorithms of 
process 25. Such formatting can and typically will 
5 include, inter alia, addition of a unique sequence 

identifier, either derived from the original accession 
number in genomic sequence database 100, or newly 
applied, and can further include additional annotation. 
Formatting can include conversion from one to another 

10 sequence listing standard, such as conversion to or 
from FAS TA or the like, depending upon the input 
expected by the subsequent process. 

Preprocessing, which can be optional 
depending upon the function desired to be identified 

15 and the informational requirements of the methods for 
effecting such identification, is followed by sequence 
processing 25, where sequences with the desired 
function are identified within the genomic sequence. 

As mentioned above, such functions can 

20 include, but are not limited to, encoding protein, 

regulating transcription, regulating message transport 
after transcription, regulating message splicing after 
transcription, regulating message degradation after 
transcription, contributing to or controlling 

25 chromosomal somatic recombination, contributing to 
chromosomal stability or movement, contributing to 
allelic exclusion or X chromosome inactivation, and the 
like . 

Where the function specified is protein 
30 coding, the above-described process of the present 
invention can be used rapidly and efficiently to 
identify individual exons in genomic sequence. 

As discussed below, and further described in 
detail in commonly owned and copending U.S. provisional 
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application nos. 60/207,456, filed May 26, 2000; 
60/234,687, filed September 21, 2000; 60/236,359, filed 
September 27, 2000; in commonly owned and copending 
U.K. patent application no. 24263.6, filed October 4, 
5 2000; and in commonly owned and copending PCT 

applications filed January 29, 2001 (attorney docket 
nos. PB 0004 WO 1, for "Human genome-derived single 
exon nucleic acid probes useful for analysis of gene 
expression in human heart"; PB 0004 WO 2, for "Human 

10 genome-derived single exon nucleic acid probes useful 
for analysis of gene expression in human brain"; PB 
0004 WO 3, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in human adult liver"; PB 0004 WO 4, for 

15 "Human genome-derived single exon nucleic acid probes 
useful for analysis of gene expression in human fetal 
liver"; PB 0004 WO 5, for "Human genome-derived single 
exon nucleic acid probes useful for analysis of gene 
expression in human lung"; PB 0004 WO 6, "Human genome- 

20 derived single exon nucleic acid probes useful for 
analysis of gene expression in human bone marrow"; 
PB 0004 WO 7, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in human placenta"; PB 0004 WO 8, for "Human 

25 genome-derived single exon nucleic acid probes useful 
for analysis of gene expression in BT 474 cells"; 
PB 0004 WO 9, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in HBL 100 cells"; PB 0004 WO 10, for "Human 

30 genome-derived single exon nucleic acid probes useful 
for analysis of gene expression in Hela cells"), the 
disclosures of which are incorporated herein by 
reference in their entirety, we have used the methods 
and apparatus of the present invention to identify more 
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than 15,000 exons in human genomic sequence whose 
expression we have confirmed in at least one human 
tissue or cell type. Fully two-thirds of the exons 
belong to genes that were not at the time of our 
5 discovery represented in existing public expression 

(EST, cDNA) databases, making the methods and apparatus 
of the present invention extremely powerful tools for 
novel gene discovery. 

And as further mentioned below and described 

10 in detail in commonly owned and copending U.S. patent 
application no. 09/632,366, filed August 3, 2000, the 
disclosure of which is incorporated herein by reference 
in its entirety, the genome-derived single exon probes 
and microarrays of the present invention prove 

15 exceedingly useful in the high throughput 

identification of a large variety of alternative splice 
events in eukaryotic cells and tissues. 

To identify such individual exons from 
genomic sequence, process 25 is used to identify 

20 putative coding regions. Two exemplary approaches 
useful in process 25 for identifying sequence that 
encodes putative genes are gene prediction and 
comparative sequence analysis. 

Gene prediction can be performed using any of 

25 a number of algorithmic methods, embodied in one or 
more software programs, that identify open reading 
frames (ORFs) using a variety of heuristics, such as 
GRAIL, DICTION, GENS CAN, and GENEFINDER. 

Comparative sequence analysis similarly can 

30 be performed using any of a variety of known programs 
that identify regions with lower sequence variability. 

An advantage of comparative sequence analysis 
is that genomic sequence can be input into process 200 
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that is less comprehensive and/or of lesser quality 
than that required by gene prediction programs. 

We have, for example, recently used 
comparative sequence analysis to identify sequences 
5 that are orthologous as between human and mouse 

genomes, and output the mouse sequences so identified 
( "similons" ) into process 300; this has permitted us to 
identify, and then to identify expression of, novel 
mouse exons and genes. As is well known in the 

10 pharmaceutical arts, genes identified in model systems 
provide targets for assessing the value of targets for 
therapeutic intervention and screening for and 
assessing agents that interact with those targets. 

As further described in Example 1, below, 

15 gene prediction software programs yield a range of 
results. For the newly accessioned human genomic 
sequence input in Example 1, for example, GRAIL 
identified the greatest percentage of genomic sequence 
as putative coding region, 2% of the data analyzed; 

20 GENEFINDER was second, calling 1%; and DICTION yielded 
the least putative coding region, with 0.8% of genomic 
sequence called as coding region. 

Increased reliability can be obtained when 
consensus is required among several such methods. 

25 Although discussed herein particularly with respect to 
exon calling, consensus among methods will in general 
increase reliability of predicting other functions as 
well . 

Thus, as indicated by query 2 6, sequence 
30 processing 25, optionally with preprocessing 24, can be 
repeated with a different method, with consensus among 
such iterations determined and reported in process 27. 

Process 27 compares the several outputs for a 
given input genomic sequence and identifies consensus 
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among the separately reported results. The consensus 
itself, as well as the sequence meeting that consensus, 
is then stored in process 29a, displayed in process 
29b, and/or output to process 300 for subsequent 
5 identification of a subset thereof suitable for assay. 
Multiple levels of consensus can be 
calculated and reported by process 27. 

For example, as further described in 
Example 1, infra, process 27 can report consensus as 

10 between all specific pairs of methods of gene 

prediction, as consensus among any one or more of the 
pairs of methods of gene prediction, or as among all of 
the gene prediction algorithms used. Thus, in Example 
1, process 27 reported that GRAIL and GENEFINDER 

15 programs agreed on 0.7% of genomic sequence, that GRAIL 
and DICTION agreed on 0.5% of genomic sequence, and 
that the three programs together agreed on 0.25% of the 
data analyzed. Put another way, 0.25% of the genomic 
sequence was identified by all three of the programs as 

2 0 containing putative coding region. 

As another example, three of the four gene 
prediction algorithms that we presently 
use - GENEFINDER, GENS CAN, and GRAIL - predict frame 
information in addition to the position of exons. If 

25 there is overlap in position and frame of the predicted 
exons, even if not complete identity, the predicted 
exons are merged in process 27 to generate the largest 
possible consensus coding region. The process is 
iterated until all possible overlaps have been merged. 

30 This approach reduces the mean number of exons present 
in each amplicon, and is preferred in generating 
exon-specif ic probes useful for detecting exon 
elongation and exon truncation alternative splice 
events . 
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Furthermore, consensus can be required among 
different approaches to identifying a chosen function. 

For example, if the function desired to be 
identified is coding of protein sequence, and a first 
5 used approach to exon calling is gene prediction, the 
process can be repeated on the same input sequence, or 
subset thereof, with another approach, such as 
comparative sequence analysis. In such a case, where 
comparative sequence analysis follows gene prediction, 

10 the comparison can be performed not only on genomic 
nucleic acid sequence, but additionally or 
alternatively can be performed on the predicted amino 
acid sequence translated from exons prior-identified by 
the gene prediction approach. 

15 Although shown as an iterative process, the 

multiple analyses required to achieve consensus can be 
done in series, in parallel, or some combination 
thereof . 

Predicted functional sequence, optionally 
20 representing a consensus among a plurality of methods 
and approaches for determination thereof, is passed to 
process 300 for identification of a subset thereof for 
functional assay. 

Where the function sought to be identified is 
25 protein coding, process 300 is used to identify a 

subset thereof suitable for experimental verification 
by physical and/or bioinf ormatic approaches. 

Where the goal is the identification and 
confirmation of expression of only a single exon of 
30 gene — for example, to provide a gene-specific 

probe — putative exons identified in process 200 can be 
classified, or binned, bioinf ormatically into putative 
genes. This binning can be based inter alia upon 
consideration of the average number of exons/gene in 
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the species chosen for analysis, upon density of exons 
that have been called on the genomic sequence, and 
other empirical rules; the putative gene structure is 
also provided by various of these gene prediction 
5 programs. Thereafter, one or more among the exons can 
be chosen for subsequent use in gene expression assay. 

Where the goal is, instead, the 
identification and confirmation of expression of all, 
or of a plurality, of the exons of a gene — as is 

10 desired for detection of alternative splice events, as 
further described in commonly owned and copending U.S. 
patent application serial no. 09/632,366, filed August 
3, 2000, the disclosure of which is incorporated herein 
by reference in its entirety — putative exons 

15 identified in process 200 can be classified, or binned, 
bioinf ormatically into putative genes. Thereafter, all 
of the exon-specif ic exons can be chosen for subsequent 
confirmation in gene expression assay. 

Where such subsequent gene expression assay 

20 uses amplified nucleic acid, considerations such as 

desired amplicon length, primer synthesis requirements, 
putative exon length, sequence GC content, existence of 
possible secondary structure, and the like can be used 
to identify and select those exons that appear most 

25 likely successfully to amplify. Where subsequent gene 
expression assay relies upon nucleic acid 
hybridization, whether or not using amplified product, 
further considerations involving hybridization 
stringency can be applied to identify that subset of 

30 sequences that will most readily permit sequence- 
specific discrimination at a chosen hybridization and 
wash stringency. One particular such consideration is 
avoidance of putative exons that span repetitive 
sequence; such sequence can hybridize spuriously to 
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nonspecific message, reducing specific signal in the 
hybridization . 

For bioinf ormatic assay, there are fewer 
constraints on the sequences that can be tested 
5 experimentally, and in this latter case therefore 
process 300 can output the entirety of the input 
sequence . 

The subset of sequences identified by 
process 300 as suitable for use in assay is then used 
10 in process 400 to create the physical and/or 

informational substrate for experimental verification 
of the predictions made in process 200, and thereafter 
to assay those substrates. 

Where the goal is to identify protein coding 
15 regions in genomic sequence, the expression of the 

sequences predicted to encode protein is verified in 
process 400 . 

Thus, in another aspect, the present 
invention provides methods and apparatus for verifying 
20 the expression of putative exons identified within 
genomic sequence. In particular, the invention 
provides methods for verifying gene expression in which 
expression of predicted exons is measured and confirmed 
using a novel type of nucleic acid microarray, the 
25 genome-derived single exon nucleic acid microarray of 
the present invention. 

According to one embodiment of this aspect, 
predicted exons are amplified from genomic DNA. 

Amplification can be performed using the 
30 polymerase chain reaction (PCR) . Although PCR is 

conveniently used, other amplification approaches, such 
as rolling circle amplification, can also be used. 

Amplification schemes can be designed to 
capture the entirety of each predicted exon in an 
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amplicon with minimal additional (that is, flanking 
intronic or intergenic) sequence. Because exons 
predicted from genomic sequence using the methods of 
the present invention differ in length, such an 
5 approach results in amplicons of varying length. 

However, we have found that most exons 
predicted from human genomic sequence are shorter than 
500 bp in length. Although amplicons of at least about 
75 base pairs, more preferably at least about 100 base 
10 pairs, even more preferably at least about 200 base 
pairs can be immobilized as probes on nucleic acid 
microarrays, our early experimental results using the 
methods of the present invention suggested that longer 
amplicons, at least about 400 base pairs, more 
15 preferably about 500 base pairs, are more effectively 
immobilized on glass slides or other prepared surfaces. 

Although we had suspected that the intronic 
and intergenic material flanking putative exons in such 
longer amplicons might cause interference with 
20 exon-specif ic hybridization during microarray 

experiments, we have found instead, to our surprise, 
that the ratio of expression of any such probe as 
between an experimental tissue (or cell type) and a 
control tissue is not significantly affected by the 
25 presence in the probes of sequence that does not 
contribute to hybridization to message or cDNA. 

Equally surprising, the art had suggested 
that single exon probes would not provide sufficient 
signal intensity for high stringency hybridization 
30 analyses. Although low stringency hybridization 

conditions have been designed that permit informative 
hybridization to highly redundant oligonucleotide-based 
microarrays, it was believed that the high stringency 
hybridization conditions typically used for EST-based 
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microarrays would not be usable with single exon 
probes. We have found, surprisingly, that single-exon 
probes provide adequate signal at high stringency. 

As a result, we have found that we are 
readily able to use genome-derived amplification 
products having a single exon flanked by intergenic 
and/or intronic sequence to confirm the expression of 
bioinf ormatically predicted exons . 

To the extent that chemical synthesis methods 
permit oligonucleotides to be generated of sufficient 
length to encompass an exon, such oligonucleotides can 
be used as probes in lieu of amplified material. At 
present, however, amplified products can be generated 
that exceed the reasonable size limit of chemically 
synthesized oligonucleotides; amplification thus more 
readily permits probes to be generated that have single 
exons flanked by intronic and/or intergenic sequence. 

Probes having flanking intergenic and/or 
intronic sequence permit a wider range of alternative 
splice events to be detected than do probes that 
contain only exonic sequence. For example, exon 
extension would be detectable with such probes as an 
increase in signal intensity: we have found a 
near-linear relationship between signal intensity and 
length of hybridizing sequence. And when used to assay 
heteronuclear, i.e., immature mRNA, probes having 
intronic and/or intergenic flanking sequence permit a 
wider variety of events to be assessed. 

Furthermore, certain advantages derive from 
application to the microarray of amplicons of defined 
size . 

Therefore, amplification schemes can 
alternatively, and preferably, be designed to amplify 
regions of defined size, preferably at least about 300 
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bp, more preferably at least about 400 bp, most 
preferably about 500 bp, centered about each predicted 
exon. Such an approach results in a population of 
amplicons of limited size diversity, but that typically 
5 contain intronic and/or intergenic nucleic acid in 
addition to, and flanking, the putative exon. 

Conversely, somewhat fewer than 10% of exons 
predicted from human genomic sequence according to the 
methods of the present invention exceed 500 bp in 

10 length. Portions of such longer exons, preferably at 
least about 300 bp, more preferably at least about 400 
bp, most preferably about 500 bp, can be amplified. 
However, in our early experiments we found that the 
percentage success at amplifying pieces of such exons 

15 is low, and that such putative exons are more 

effectively amplified when larger fragments, at least 
about 1000 bp, typically at least about 1500 bp, and 
even as large as 2000 bp are amplified. Further 
routine optimization of the PCR reaction would permit 

20 500 bp portions of the longer exons to be amplified. 

For amplification, the putative exons 
selected in process 300 are input into one or more 
primer design programs, such as PRIMER3 (available 
online for use at 

25 http://www-genome.wi.mit.edu/cgi-bin/primer/ ), with a 
goal of amplifying at least about 500 base pairs of 
genomic sequence centered within or about exons 
predicted to be no more than about 500 bp, or at least 
about 1000 - 1500 bp of genomic sequence for exons 

30 predicted to exceed 500 bp in length, and the primers 
synthesized by standard techniques. Primers with the 
requisite sequences can be purchased commercially or 
synthesized by standard techniques. 



Conveniently, a first predetermined sequence 
can be added commonly to each exon-specif ic 5 T primer 
and a second, typically different, predetermined 
sequence commonly added to each 3 1 exon-unique primer. 
5 This serves to immortalize the amplicon: that is, it 
serves to permit further amplification of any amplicon 
using a single set of primers complementary 
respectively to the common 5 T and common 3 1 sequence 
elements. The presence of these "universal" priming 

10 sequences further facilitates later sequence 

verification, providing a sequence common to all 
amplicons at which to prime sequencing reactions. The 
common 5 T and 3 T sequences can further serve to add a 
cloning site should any of the exons warrant further 

15 study. 

Such predetermined sequence is usefully at 
least about 10 nt in length, typically at least about 
12 nt, more typically about 15 nt in length, and 
usually does not exceed about 25 nt in length. The 

20 "universal" priming sequences used in the examples 

presented infra were each 16 nt long, and are further 
described in commonly owned and copending U.S. patent 
application serial no. 09/608,408, filed June 30, 2000, 
the disclosure of which is incorporated herein by 

25 reference in its entirety. 

The genomic DNA to be used as substrate for 
amplification will come from the eukaryotic species 
from which the genomic sequence data had originally 
been obtained, or a closely related species, and can 

30 conveniently be prepared by well known techniques from 
somatic or germline tissue or cultured cells of the 
organism. See, e.g., Short Protocols in Molecular 
Biology : A Compendium of Methods from Current 
Protocols in Molecular Biology , Ausubel et al . (eds.), 



4 th edition (April 1999), John Wiley & Sons (ISBN: 
047132938X) and Maniatis et al., Molecular Cloning : A 
Laboratory Manual , 2 nd edition (December 1989), Cold 
Spring Harbor Laboratory Press (ISBN: 0879693096), the 
disclosures of which are incorporated herein by 
reference in their entireties. Many such prepared 
genomic DNAs are available commercially, with the human 
genomic DNAs additionally having certification of donor 
informed consent. 

After partial purification, as by size 
exclusion spin column or adsorption to glass, with or 
without confirmation as to amplicon quality as by gel 
electrophoresis, each amplicon (single exon probe) is 
disposed in an array upon a support substrate. 

Methods for creating microarrays by 
deposition and fixation of nucleic acids onto support 
substrates are well known in the art. Reviewed in 
Schena (ed.), DNA Microarravs : A Practical Approach 
(Practi cal Approach Series ) , Oxford University Press 
(1999) (ISBN: 0199637768); Nature Genet. 
21 (1) (suppl) : 1 - 60 (1999); Schena (ed.), Microarrav 
Biochip: Tools and Technology , Eaton Publishing 
Company/BioTechniques Books Division (2000) (ISBN: 
1881299376), the disclosures of which are incorporated 
herein by reference in their entireties. 

Typically, the support substrate can be 
glass, although other materials, such as amorphous 
silicon, crystalline silicon, or plastics, can be used. 
Such plastics include polymethylacrylic, polyethylene, 
polypropylene, polyacrylate, polymethylmethacrylate, 
polyvinylchloride, polytetraf luoroethylene, 
polystyrene, polycarbonate, polyacetal, polysulfone, 
celluloseacetate, cellulosenitrate, nitrocellulose, or 
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mixtures thereof. Typically, the support can be 
rectangular, although other shapes, particularly 
circular disks and even spheres, present certain 
advantages. Particularly advantageous alternatives to 
5 glass slides as support substrates for array of nucleic 
acids are optical discs, as described in Demers, 
"Spatially Addressable Combinatorial Chemical Arrays in 
CD-ROM Format, " international patent publication 
WO 98/12559, incorporated herein by reference in its 
10 entirety. 

The amplified nucleic acids can be attached 
covalently to a surface of the support substrate or, 
more typically, applied to a derivatized surface in a 
chaotropic agent that facilitates denaturation and 

15 adherence by presumed noncovalent interactions, or some 
combination thereof . 

Robotic spotting devices useful for arraying 
nucleic acids on support substrates can be constructed 
using public domain specifications (The MGuide, version 

20 2.0, http : //cmgm. Stanford. edu/pbrown/mguide/ 

index.html), or can conveniently be purchased from 
commercial sources (MicroArray Genii Spotter and 
MicroArray Genlll Spotter, Molecular Dynamics, Inc., 
Sunnyvale, CA) . Spotting can also be effected by 

25 printing methods, including those using ink jet 
technology. 

As is well known in the art, microarrays 
typically also contain immobilized control nucleic 
acids. For controls useful in providing measurements 

30 of background signal for the genome-derived single exon 
microarrays of the present invention, a plurality of 
E. coll genes can readily be used. As further 
described in Example 1, 16 or 32 E. coli genes suffice 
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to provide a robust measure of nonspecific 
hybridization in such microarrays. 

As is well known in the art, the amplified 
product disposed in arrays on a support substrate to 
5 create a nucleic acid microarray can consist entirely 
of natural nucleotides linked by phosphodiester bonds, 
or alternatively can include either nonnative 
nucleotides, alternative internucleotide linkages, or 
both, so long as complementary binding can be obtained 

10 in the hybridization reaction. If enzymatic 

amplification is used to produce the immobilized 
probes, the amplifying enzyme will impose certain 
further constraints upon the types of nucleic acid 
analogs that can be generated, 

15 Although particularly described herein as 

using high density microarrays constructed on planar 
substrates, the methods of the present invention for 
confirming the expression of exons predicted from 
genomic sequence can use any of the known types of 

20 microarrays as herein defined, including microarrays on 
nonplanar, nonunitary, distributed substrates, such as 
the nonplanar, bead-based microarrays as are described 
in Brenner et al . , Proc. Natl. Acad. Scl. USA 
97 (4) : 166501670 (2000); U.S. Patent No. 6, 057, 107; and 

25 U.S. Patent No. 5,736,330, the disclosures of which are 
incorporated herein by reference in their entireties. 
In theory, a packed collection of such beads provides 
in aggregate a higher density of nucleic acid probe 
than can be achieved with spotting or lithography 

30 techniques on a single planar substrate. 

In addition, gene expression can be confirmed 
using hybridization to lower density arrays, such as 
those constructed on membranes, such as nitrocellulose, 
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nylon, and positively-charged derivatized nylon 
membranes . 

Planar microarrays on solid substrates, 
however, provide certain useful advantages, including 
5 compatibility with existing readers. For example, each 
standard microscope slide can include at least 1000, 
typically at least 2000, preferably 5000 or more, and 
up to 19,000 or more nucleic acid probes of discrete 
sequence . 

10 Each putative gene can be represented in the 

array by a single predicted exon or by a plurality of 
exons predicted to belong to the same gene. And as is 
well known in the art, each probe of defined sequence, 
representing a single predicted exon, can be deposited 

15 in a plurality of locations on a single microarray to 
provide redundancy of signal. 

The genome-derived single exon microarrays 
described above are an important aspect of the present 
invention, and differ in several fundamental and 

20 advantageous ways from microarrays presently used in 

the gene expression art, including (1) those created by 
deposition of mRNA-derived nucleic acids, (2) those 
created by in situ synthesis of oligonucleotide probes, 
and (3) those constructed from yeast genomic DNA. 

25 Most nucleic acid microarrays that are in use 

for study of eukaryotic gene expression have as 
immobilized probes nucleic acids that are derived — 
either directly or indirectly — from expressed message. 
It is common, for example, for such microarrays to be 

30 derived from cDNA/EST libraries, either from those 

previously described in the literature, such as those 
from the I.M.A.G.E. consortium, Lennon et al . , "The 
I.M.A.G.E. Consortium: an Integrated Molecular Analysis 
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of Genomes and Their Expression, Genomics 33(1): 151-2 
(1996), or from the de novo construction of "problem 
specific" libraries targeted at a particular biological 
question, R.S. Thomas et al . , Toxicologist 54:68-69 
5 (2000), incorporated herein by reference in their 

entireties. Such microarrays are herein collectively 
denominated "EST microarrays". 

Such EST microarrays by definition can 
measure expression only of those genes found in EST 

10 libraries, which we show herein (see infra) to 

represent only a fraction of expressed genes. Thus, as 
further discussed in Example 1, infra, fully 2/3 of 
genes identified from newly-accessioned human genomic 
sequence data by the methods of the present 

15 invention — for which expression was subsequently 
confirmed using the methods and apparatus of the 
present invention — do not appear in EST or other 
expression databases, and could not, therefore, have 
been represented as probes on an EST microarray. 

20 Furthermore, EST and cDNA libraries — and 

thus microarrays based thereupon — are biased by the 
tissue or cell type of message origin. 

In addition, representation of a message in 
an EST and/or cDNA library depends upon the successful 

25 reverse transcription, optionally but typically with 
subsequent successful cloning, of the message. This 
introduces substantial bias into the population of 
probes available for arraying in EST microarrays. For 
example, as we show in the examples, infra, the subset 

30 of genes identified from genomic sequence by the 

methods of the present invention that had previously 
been accessioned in EST or other expression 
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databases are biased toward genes with higher 
expression levels . 

In contrast, neither reverse transcription 
nor cloning is required to produce the probes arrayed 
5 on the genome-derived single exon microarrays of the 
present invention. And although the ultimate 
deposition of a probe on the genome-derived single exon 
microarray of the present invention depends upon a 
successful amplification from genomic material, a 

10 priori knowledge of the sequence of the desired 

amplicon affords greater opportunity to recover any 
given probe sequence recalcitrant to amplification than 
is afforded by the requirement for successful reverse 
transcription and cloning of unknown message in EST 

15 approaches. Furthermore, if the sequence cannot be 
amplified, the sequence can at times be chemically 
synthesized in its entirety for use in the present 
invention. 

Thus, the genome-derived single exon 
20 microarrays of the present invention present a far 
greater diversity of probes for measuring gene 
expression, with far less bias, than do EST microarrays 
presently used in the art. 

As a further consequence of their ultimate 
25 origin from expressed message, the probes in EST 
microarrays often contain poly-A (or complementary 
poly-T) stretches derived from the poly-A tail of 
mature mRNA. These homopolymeric stretches contribute 
to cross-hybridization, that is, to a spurious signal 
30 occasioned by hybridization to the homopolymeric tail 
of a labeled cDNA that lacks sequence homology to the 
gene-specific portion of the probe. 

In contrast, the probes arrayed in the 
genome-derived single exon microarrays of the present 
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invention lack homopolymeric stretches derived from 
message polyadenylation, and thus can provide more 
specific signal. Typically, at least about 50% of the 
probes on the genome-derived single exon microarrays of 
5 the present invention lack homopolymeric regions 

consisting of A or T, where a homopolymeric region is 
defined for purposes herein as stretches of 25 or more, 
typically 30 or more, identical nucleotides. More 
typically, at least about 60%, even more typically at 

10 least about 75%, of probes on the genome-derived single 
exon microarrays of the present invention lack such 
homopolymeric stretches ♦ 

A further distinction, which also affects the 
specificity of hybridization, is occasioned by the 

15 typical derivation of EST microarray probes from cloned 
material. Because much of the probe material disposed 
as probes on EST microarrays is excised or amplified 
from plasmid, phage, or phagemid vectors, EST 
microarrays typically include a fair amount of vector 

20 sequence, more so when the probes are amplified, rather 
than excised, from the vector. 

In contrast, the vast majority of probes in 
the genome-derived single exon microarrays of the 
present invention contain no prokaryotic or 

25 bacteriophage vector sequence, having been amplified 
directly or indirectly from genomic DNA. Typically, 
therefore, at least about 50%, more typically at least 
about 60%, 70%, and even 80% or more of individual 
exon-including probes disposed on a genome-derived 

30 single exon microarray of the present invention lack 

vector sequence, and particularly lack sequences drawn 
from plasmids and bacteriophage. Preferably, at least 
about 85%, more preferably at least about 90%, most 
preferably more than 90% of exon-including probes in 
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the genome-derived single exon microarray of the 
present invention lack vector sequence. With attention 
to removal of vector sequences through preprocessing 
24, percentages of vector-free exon-including probes 
5 can be as high as 95 - 99%. The substantial absence of 
vector sequence from the genome-derived single exon 
microarrays of the present invention results in greater 
specificity during hybridization, since spurious cross- 
hybridization to a probe vector sequence is reduced. 

10 As a further consequence of excision or 

amplification of probes from vectors in construction of 
EST microarrays, the probes arrayed thereon often 
contain artificial sequence, derived from vector 
polylinker multiple cloning sites, at both 5 1 and 3 1 

15 ends. The probes disposed upon the genome-derived 
single exon microarrays need have no such artificial 
sequence appended thereto. 

As mentioned above, however, the exon- 
specific primers used to amplify putative exons can 

20 include artificial sequences, typically 5 T to the exon- 
specific primer sequence, useful for "universal" (that 
is, independent of exon sequence) priming of subsequent 
amplification or sequencing reactions. When such 
"universal" 5 T and/or 3 T priming sequences are appended 

25 to the amplification primers, the probes disposed upon 
the genome-derived single exon microarray will include 
artificial sequence similar to that found in EST 
microarrays. However, the genome-derived single exon 
microarray of the present invention can be made without 

30 such sequences, and if so constructed, presents an even 
smaller amount of nonspecific sequence that would 
contribute to nonspecific hybridization. 

Yet another consequence of typical use of 
cloned material as probes in EST microarrays is that 



such microarrays contain probes that result from 
cloning artifacts, such as chimeric molecules 
containing coding region of two separate genes. 
Derived from genomic material, typically not thereafter 
cloned, the probes of the genome-derived single exon 
microarrays of the present invention lack such cloning 
artifacts, and thus provide greater specificity of 
signal in gene expression measurements, 

A further consequence of the cloned origin of 
probes on many EST microarrays is that the individual 
probes often have disparate sizes, which can cause the 
optimal hybridization stringency to vary among probes 
on a single microarray. In contrast, as discussed 
above, the probes arrayed on the genome-derived single 
exon microarrays of the present invention can readily 
be designed to have a narrow distribution in sizes, 
with the range of probe sizes no greater than about 10% 
of the average size, typically no greater than about 5% 
of the average probe size. 

Because of their origin from fully- or 
partially-spliced message, probes disposed upon EST 
arrays will often include multiple exons . The 
percentage of such exon-spanning probes in an EST 
microarray can be calculated, on average, based upon 
the predicted number of exons/gene for the given 
species and the average length of the immobilized 
probes. For human genes, the near-complete sequence of 
human chromosome 22, Dunham et al . , Nature 
402 (6761) : 489-95 (1999), predicts that human genes 
average 5.5 exons/gene. Even with probes of 200 - 
500 bp, the vast majority of human EST microarray 
probes include more than one exon. 

In contrast, by virtue of their origin from 
algorithmically identified exons in genomic sequence, 
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the probes in the genome-derived single exon 
microarrays of the present invention can comprise 
individual exons, which provides the ability, as 
further discussed in commonly owned and copending U.S. 
5 patent application serial no. 09/632,366, filed 

August 3, 2000, incorporated herein by reference in its 
entirety, to detect and to characterize the expression 
of splice variants. 

Although the presence of multiexon probes 

10 will not interfere with the ability to confirm 

expression of predicted exons in a first level screen, 
it is preferred that at least about 50%, typically at 
least about 60%, even more typically at least about 70% 
of probes disposed on the genome-derived microarray of 

15 the present invention consist of, or include, no more 
than one exon. In preferred embodiments, at least 
about 75%, more preferably at least about 80%, 85%, 
90%, 95%, and even 99% of probes in the genome-derived 
microarrays of the present invention consist of, or 

20 include, no more than one exon. 

Although, in the most preferred embodiments, 
at least about 95%, and even at least about 99% of 
probes in the genome-derived microarray consist of, or 
include, no more than one exon, we have found that our 

25 early bioinf ormatic parameters typically produce, at 
this stage of analysis, about 10% of probes that 
potentially contain two exons. We expect that some 
fraction of these probes will prove to encode only a 
single exon, and that further optimization of our 

30 bioinformatic approach will reduce the percentage of 
probes having more than one potential exon. 

Further distinguishing the genome-derived 
single exon microarrays of the present invention from 
the EST arrays in the art, the exons that are 
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represented in EST microarrays are often biased toward 
the 3' or 5 T end of their respective genes, since 
sequencing strategies used for EST identification are 
so biased. In contrast, no such 3 1 or 5 f bias 
5 necessarily inheres in the selection of exons for 
disposition on the genome-derived single exon 
microarrays of the present invention. 

Conversely, the probes provided on the 
genome-derived single exon microarrays of the present 

10 invention typically, but need not necessarily, include 
intronic and/or intergenic sequence that is absent from 
EST microarrays, which are derived from mature mRNA. 
As above mentioned, such inclusion, although not 
mandatory, is advantageous, particularly in use of the 

15 probes for detection of alternative splice events. 
Typically, therefore, at least about 50%, more 
typically at least about 60%, and even more typically 
at least about 70% of the exon-including probes on the 
genome-derived single exon microarrays of the present 

2 0 invention include sequence drawn from noncoding 

regions. In some embodiments, at least about 80 % , mo r e 
typically at least about 85%, 90%, 91%, 92%, 93%, 94%, 
95%, 96%, 97%, 98%, and even 99% or more of exon- 
including probes on the genome-derived single exon 

25 microarrays of the present invention will include 
sequence drawn from noncoding regions. 

The genome-derived single exon microarrays of 
the present invention are also quite different from in 
situ synthesis microarrays, where probe size is 

30 severely constrained by limitations of the 

photolithographic or other in situ synthesis processes. 

Typically, probes arrayed on in situ 
synthesis microarrays are limited to a maximum of about 



- 46 - 



25 bp. As a well known consequence, hybridization to 
such chips must be performed at low stringency. In 
order, therefore, to achieve unambiguous sequence- 
specific hybridization results, the in situ synthesis 
5 microarray requires substantial redundancy, with 

concomitant programmed arraying for each probe of probe 
analogues with altered (i.e., mismatched) sequence. 

In contrast, the longer probe length of the 
genome-derived single exon microarrays of the present 

10 invention allows much higher stringency hybridization 
and wash. Typically, therefore, exon-including probes 
on the genome-derived single exon microarrays of the 
present invention average at least about 100 bp, more 
typically at least about 200 bp, preferably at least 

15 about 250 bp, even more preferably about 300 bp, 400 
bp, or in preferred embodiments, at least about 500 bp 
in length. By obviating the need for substantial probe 
redundancy, this approach permits a higher density of 
probes for discrete exons or genes to be arrayed on the 

20 microarrays of the present invention than can be 
achieved for in situ synthesis microarrays. 

A further distinction is that the probes in 
in situ synthesis microarrays typically are covalently 
linked to the substrate surface. In contrast, the 

25 probes disposed on the genome-derived microarray of the 
present invention typically are, but need not 
necessarily be, bound noncovalently to the substrate. 

Furthermore, the short probe size on in situ 
microarrays causes large percentage differences in the 

30 melting temperature of probes hybridized to their 

complementary target sequence, and thus causes large 
percentage differences in the theoretically optimum 
stringency across the array as a whole. 
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In contrast, the larger probe size in the 
microarrays of the present invention create lower 
percentage differences in melting temperature across 
the range of arrayed probes. 
5 A further significant advantage of the 

microarrays of the present invention over in situ 
synthesized arrays is that the quality of each 
individual probe can be confirmed before deposition. 
In contrast, the quality of probes cannot be assessed 

10 on a probe-by-probe basis for the in situ synthesized 
microarrays presently being used. 

The genome-derived single exon microarrays of 
the present invention are also distinguished over, and 
present substantial benefits over, the genome-derived 

15 microarrays from lower eukaryotes such as yeast. See, 
e.g., Lashkari et al . , Proc. Natl. Acad. Sci. USA 
94 : 13057-13062 (1997) . 

Only about 220 - 250 of the 6100 or so 
nuclear genes in Saccharomyces cerevisiae — that is, 

20 only about 4 to 5% — have standard, spliceosomal , 
introns, Lopez et al . , Nucl. Acids Res. 28:85-86 
(2000); Spingola et al . , RNA 5(2):221-34 (1999), 
permitting the ready amplification and disposition of 
single-exon amplicons on such microarray without the 

25 requirement for antecedent use of gene prediction 
and/or comparative sequence analyses. 

A significant aspect of the present invention 
is the ability to identify and to confirm expression of 
predicted coding regions in genomic sequence drawn from 

30 eukaryotic organisms that have a higher percentage of 
genes having introns than do yeast such as 
Saccharomyces cerevisiae, particularly in genomic 
sequence drawn from eukaryotes in which at least about 
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10%, typically at least about 20%, more typically at 
least about 50% of protein-encoding genes have introns. 
In preferred embodiments, the methods and apparatus of 
the present invention are used to identify and confirm 
5 expression of exons of novel genes from genomic 

sequence of eukaryotes in which the average number of 
introns per gene is at least about one, more typically 
at least about two, even more typically at least about 
three or more. 

10 After the physical substrate is prepared, 

experimental verification of predicted function is 
performed. 

In a preferred embodiment of the present 
invention, where the function sought to be identified 

15 in genomic sequence is protein coding, experimental 
verification is performed by measuring expression of 
the putative exons, typically through nucleic acid 
hybridization experiments, and in particularly 
preferred embodiments, through hybridization to genome- 

20 derived single exon microarrays prepared as above 
described. 

Expression is conveniently measured and 
reported for each probe in the microarray both as a 
signal intensity and as a ratio of the expression 

25 measured relative to a control, according to techniques 
well known in the microarray art, reviewed in Schena 
(ed.), DNA Microarrays : A Practical Approach 
(Practical Approach Series ) , Oxford University Press 
(1999) (ISBN: 0199637768); Nature Genet. 

30 21 (1) (suppl) : 1 - 60 (1999); Schena (ed.), Microarray 
Biochip: Tools and Technology , Eaton Publishing 
Company/BioTechniques Books Division (2000) (ISBN: 
1881299376), the disclosures of which are incorporated 
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herein by reference in their entireties. See also 
Example 2, infra. The mRNA source for the reference 
(control) used to calculate expression ratios can be 
heterogeneous, as from a pool of multiple tissues 
5 and/or cell types or, alternatively, can be drawn from 
a homogeneous mRNA source, such as a single cultured 
cell-type . 

In Examples 1 and 2, infra, we used a pool of 
10 tissues/cell types as control. We have since 

10 observed that almost every probe that demonstrates 

expression in the control pool can readily be shown to 
be expressed in HeLa cells. Since use of a pooled 
control might mask subtle alternative splice events, we 
have used HeLa as the source of control message in more 

15 recent experiments. 

mRNA can be prepared by standard techniques, 
Short Protocols in Molecular Biology : A Compendium of 
Methods from Current Protocols in Molecular Biology , 
Ausubel et al. (eds.), 4 th edition (April 1999), John 

20 Wiley & Sons (ISBN: 047132938X) and Maniatis et al . , 
Molecular Cloning : A Laboratory Manual , 2 nd edition 
(December 1989) , Cold Spring Harbor Laboratory Press 
(ISBN: 0879693096), the disclosures of which are 
incorporated herein by reference in their entireties, 

25 or purchased commercially. The mRNA is then typically 
reverse-transcribed in the presence of labeled 
nucleotides: the index source (that in which expression 
is desired to be measured) is reverse transcribed in 
the presence of nucleotides labeled with a first label, 

30 typically a fluorophore (equivalently denominated 

f luorochrome; fluor; fluorescent dye) ; the reference 
source is reverse transcribed in the presence of a 
second label, typically a fluorophore, typically 
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f luorometrically-distinguishable from the first label. 
As further described in Example 2, infra, Cy3 and Cy5 
dyes prove particularly useful in these methods. After 
partial purification of the index and reference 
5 targets, hybridization to the probe array is conducted 
according to standard techniques, typically under a 
coverslip or in an automatic slide processing unit. 

After wash, microarrays are conveniently 
scanned using a commercial microarray scanning device, 

10 such as a Gen3 or Avalanche Scanner (Molecular 

Dynamics, Sunnyvale, CA) . Data on expression is then 
passed, with or without interim storage, to process 
500, where the results for each probe are related to 
the original sequence. 

15 Often, hybridization of target material to 

the genome-derived single exon microarray will identify 
certain of the probes thereon as of particular 
interest. Thus, it is often desirable that the user be 
able readily to obtain sufficient quantities of an 

20 individual probe, either for subsequent arrayed 

deposition upon an additional support substrate, often 
as part of a microarray having a plurality of probes so 
identified, or alternatively or additionally as a 
solitary solid-phase or solution-phase probe for 

25 further use. 

Thus, in another aspect, the present 
invention provides compositions and kits for the ready 
production of nucleic acids identical in sequence to, 
or substantially identical in sequence to, probes on 

30 the genome-derived single exon microarrays of the 
present invention . 

In one embodiment, the invention provides 
individual single exon probes in the form of 
substantially isolated and purified nucleic acid. In 
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one such embodiment the probe is provided in quantity 
sufficient to perform a hybridization reaction. 

When provided in quantity sufficient to 
perform a hybridization reaction, the probe can be in 
5 any form directly hybridizable to the target that 

contains the probe T s exon (or its complement), such as 
double stranded DNA, single-stranded DNA complementary 
to the target, single-stranded RNA complementary to the 
target, or chimeric DNA/ RNA molecules so hybridizable. 

10 The nucleic acid can alternatively or 

additionally include either nonnative nucleotides, 
alternative internucleotide linkages, or both, so long 
as complementary binding can be obtained. For example, 
probes can include phosphorothioates, 

15 methylphosphonates, morpholino analogs, and peptide 

nucleic acids (PNA) , as are described, inter alia, in 
U.S. Patent Nos. 5,142,047; 5,235,033; 5,166,315; 
5,217,866; 5,184,444; 5,861,250; international patent 
applications nos. WO 93/25706; and in Science 254:1497 

20 (1991); J. Am. Chem. Soc. 114:9677 (1992); J. Am. Chem. 
Soc. 144:1895 (1992); J. Chem. Soc. Chem. Comm. 800 
(1993); Proc. Nat. Acad. Sci. USA 90:1667 (1993); 
Intercept Ltd. 325 (1992); J. Am. Chem. Soc. 114:9677 

(1992) ; Nucleic Acids Res. 21:197 (1993); J. Chem. Soc. 
25 Chem. Commun. 518 (1993); Anti-Cancer Drug Design 8:53 

(1993) ; Nucleic Acids Res. 21:2103 (1993); Org. Proc. 
Prep. 25:457 (1993); CRC Press 363 (1992); J. Chem. 
Soc. Chem. Commun. 9:800 (1993); J. Am. Chem. Soc. 
115:6477 (1993); Nature 365:566 (1993); WO 92/20702; 

30 and WO 92/20703, the disclosures of which are 
incorporated herein by reference. 

Usefully, however, such probes are instead 
provided in a form and quantity suitable for 
amplification, such as by PCR. Although PCR is 
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conveniently used, other amplification approaches can 
be used as well, such as rolling circle amplification, 
as is described, Inter alia, in U.S. Patent Nos. 
5,854,033 and 5,714,320 and international patent 
5 publications WO 97/19193 and WO 00/15779, the 
disclosures of which are incorporated herein by 
reference in their entireties. As is well understood, 
where the probes are to be provided in a form suitable 
for amplification, the range of nucleic acid analogues 
10 and/or internucleotide linkages will be constrained by 
the requirements and nature of the amplification 
enzyme . 

Where the probe is to be provided in form 
suitable for amplification, the quantity need not be 

15 sufficient for direct hybridization for gene expression 
analysis, and need be sufficient only to function as an 
amplification template, typically at least about 1 pg, 
more typically at least about 10 pg, and usually at 
least about 100 pg or more. 

20 Each discrete amplifiable probe can also be 

packaged with amplification primers, either in a single 
composition that comprises probe template and primers, 
or in a kit that comprises such primers separately 
packaged therefrom. As above mentioned, the 

25 exon-specif ic 5 1 primers used for genomic amplification 
can have a first common sequence added thereto, and the 
exon-specif ic 3 f primers used for genomic amplification 
can have a second, different, common sequence added 
thereto, thus permitting, in this embodiment, the use 

30 of a single set of 5 T and 3 T primers to amplify any one 
of the probes. The probe composition and/or kit can 
also include buffers, enzyme, etc., required to effect 
amplification. 
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In another embodiment , only amplification 
primers are provided. The primers are sufficient to 
permit generation of the single exon probe by 
amplification from genomic DNA, which can be provided 
5 by the user. 

As mentioned above, when intended for use on 
a genome-derived single exon microarray of the present 
invention, the genome-derived single exon probes of the 
present invention will typically average at least about 

10 75 - 100 bp, more typically at least about 200 bp, 

preferably at least about 250 bp, even more preferably 
about 300 bp, 400 bp, or in preferred embodiments, at 
least about 500 bp in length, including (and typically, 
but not necessarily centered about) the exon. 

15 Furthermore, when intended for use on a genome-derived 
single exon microarray of the present invention, the 
genome-derived single exon probes of the present 
invention will typically not contain a detectable 
label . 

20 When intended for use in solution phase 

hybridization, however — that is, for use in a 
hybridization reaction in which the probe is not first 
bound to a support substrate (although the target may 
indeed be so bound) — length constraints that are 

25 imposed in microarray-based hybridization approaches 
will be relaxed, and such probes will typically be 
labeled. 

In such case, the only functional constraint 
that dictates the minimum size of such probe is that 
30 each such probe must be capable of specifically 

identifying in a hybridization reaction the exon from 
which it is drawn. In theory, a probe of as little as 
17 nucleotides is capable of uniquely identifying its 
cognate sequence in the human genome. For 



hybridization to expressed message — a subset of target 
sequence that is much reduced in complexity as compared 
to genomic sequence - even fewer nucleotides are 
required for specificity. 

Therefore, the probes of the present 
invention can include as few as 20 bp of exon, 
typically at least about 25 bp of exon, more typically 
at least about 50 bp or exon, or more. The minimum 
amount of exon required to be included in the probe of 
the present invention in order to provide specific 
signal in either solution phase or microarray-based 
hybridizations can readily be determined by routine 
experimentation using standard high stringency 
conditions . 

Such high stringency conditions are 
described, inter alia, in Short Protocols in Molecular 
Biology : A C ompendium of Methods from Current 
Protocols in Mol ecular Biology . Ausubel et al. (eds.), 
4th edition (April 1999), John Wiley & Sons (ISBN: 
047132938X) and Maniatis et al . , Molecular Cloning ! A 
Laboratory Manual, 2nd edition (December 1989) , Cold 
Spring Harbor Laboratory Press (ISBN: 0879693096), the 
disclosures of which are incorporated herein by 
reference in their entireties. 

For microarray-based hybridization, standard 
high stringency conditions can usefully be 50% 
formamide, 5X SSC, 0.2 ]iq/]il poly(dA), 0.2 ug/ul human 
cotl DNA, and 0.5 % SDS, in a humid oven at 42°C 
overnight, followed by successive washes of the 
microarray in IX SSC, 0.2% SDS at 55°C for 5 minutes, 
and then 0 . IX SSC, 0.2% SDS, at 55°C for 20 minutes. 



- 55 - 



For solution phase hybridization, standard 
high stringency conditions can usefully be aqueous 
hybridization at 65°C in 6X SSC. 

Lower stringency conditions, suitable for 
5 cross-hybridization to mRNA encoding structurally- and 
functionally-related proteins, can usefully be the same 
as the high stringency conditions but with reduction in 
temperature for hybridization and washing to room 
temperature (approximately 25°C) . 
10 When intended for use in solution phase 

hybridization, the maximum size of the single exon 
probes of the present invention is dictated by the 
proximity of other exons in genomic DNA: although each 
single exon probe can include intergenic and/or 
15 intronic material contiguous to the exon in the human 
genome, each probe of the present invention will 
typically include portions of only one exon. 

Thus, each single exon probe will include no 
more than about 25 kb of contiguous genomic sequence, 
2 0 more typically no more than about 20 kb of contiguous 
genomic sequence, more usually no more than about 
15 kb, even more usually no more than about 10 kb . 
Usually, probes that are maximally about 5 kb will be 
used, more typically no more than about 3 kb. 
25 xt will be appreciated that single stranded 

probes must be complementary in sequence to the target; 
it is well within the skill in the art to determine 
such complementary sequence and the need therefor. it 
will further be understood that double stranded probes 
30 can be used in both solution-phase hybridization and 
microarray-based hybridization if suitably denatured. 
Thus, it is an aspect of the present invention to 
provide single-stranded nucleic acid probes that have 
sequence complementary to those described herein above 



and below, and double-stranded probes one strand of 
which has sequence complementary to the probes 
described herein. 

As mentioned above, the probes can, but need 
not, contain intergenic and/or intronic material that 
flanks the exon, on one or both sides, in the same 
linear relationship to the exon that the intergenic 
and/or intronic material bears to the exon in genomic 
DNA. The probes typically do not, however, contain 
nucleic acid derived from more than one expressed exon. 

And when intended for use in solution 
hybridization, the probes of the present invention can 
usefully have detectable labels. Nucleic acid labels 
are well known in the art, and include, inter alia, 
radioactive labels, such as 3 H, 32 P, 33 p, 35 s, 125 I, 131 i; 
fluorescent labels, such as Cy3, Cy5, Cy5.5, Cy7, SYBR® 
Green and other labels described in Haugland, Handbook 
of Fluorescent Probes and Research Chemicals , 7th ed., 
Molecular Probes Inc., Eugene, OR (2000), or 
fluorescence resonance energy transfer tandem 
conjugates thereof; labels suitable for 
chemiluminescent and/or enhanced chemiluminescent 
detection; labels suitable for ESR and NMR detection; 
quantum dots; and labels that include one member of a 
specific binding pair, such as biotin, digoxigenin, or 
the like. 

The probes, either in quantity sufficient for 
hybridization or sufficient for amplification, can be 
provided in individual vials or containers, and can be 
provided dry (e.g., lyophilized) , or solvated. If 
solvated, the solution can usefully include buffers and 
salts as desired for hybridization and/or 
amplification. Furthermore, if desired to be spotted 
on a microarray, the probes can usefully be provided in 



a solution of chaotropic agent to facilitate adherence 
to the microarray support substrate. 

Alternatively, such probes can usefully be 
packaged as a plurality of such individual 
genome-derived single exon probes. 

In one embodiment of this aspect, a small 
quantity of each probe is disposed, typically without 
attachment to substrate, in a spatially-addressable 
ordered set, typically one per well of a microtiter 
dish. Although a 96 well microtiter plate can be used, 
greater efficiency is obtained using higher density 
arrays, such as are provided by microtiter plates 
having 384, 864, 1536, 3456, 6144, or 9600 wells. And 
although microtiter plates having physical depressions 
(wells) are conveniently used, any device that permits 
addressable withdrawal of reagent from fluidly- 
noncommunicating areas can be used. 

Each of the probes of the ordered set can be 
provided in any of the forms that are described above 
with respect to the probes as individually packaged. 

As above mentioned, the exon-specif ic 
5 T primers used for genomic amplification can have a 
first common sequence added thereto, and the exon- 
specif ic 3 T primers used for genomic amplification can 
have a second, different, common sequence added 
thereto, thus permitting, in certain embodiments, the 
use of a single set of 5 T and 3 f primers to amplify any 
one of the probes from the amplifiable ordered set. 

Such collections of genome-derived single 
exon probes can usefully include a plurality of probes 
chosen for a common attribute, such as common 
expression in a given tissue, cell type, developmental 
stage, disease state, or the like. 
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In such defined subsets, typically at least 
50% of the probes will have the common attribute, such 
as expression in the defined tissue or cell type. More 
typically, at least about 60% of the probes will be 
5 expressed in the defined tissue, even more typically at 
least about 75%, and preferably at least about 80%, 
85%, or, in preferred embodiments, at least about 90%, 
and even 95% or more of the probes will have the common 
attribute, such as expression in the defined tissue or 
10 cell type. 

Analogously, the invention provides, in 
another aspect, genome-derived single-exon nucleic acid 
microarrays having a plurality of probes chosen for a 
common attribute, such as common expression in a given 

15 tissue, cell type, developmental stage, disease state, 
or the like. 

These "subset-defined" genome-derived single 
exon microarrays can be distinguished from the "first 
iteration" genome-derived single exon microarrays of 

20 the present invention, i.e., from those that are used 
to confirm expression of predicted exons, by the 
percentage of probes that are known to have a common 
attribute, such as expression in a defined tissue or 
cell type. On such "subset-defined" microarrays, 

25 typically at least 50% of the probes will have the 

common attribute, typically expression in the defined 
tissue or cell type. More typically, at least about 
60% of the probes will be expressed in the defined 
tissue, even more typically at least about 75%, and 

30 preferably at least about 80%, 85%, or, in preferred 
embodiments, at least about 90%, and even 95% or more 
of the probes will have the common attribute, such as 
expression in the defined tissue or cell type. 



When used for gene expression analysis, the 
"defined subset" genome-derived single exon microarrays 
provide greater physical informational density than do 
the genome-derived single exon microarrays that have 
lower percentages of probes known to be expressed 
commonly in the tested tissue. At a fixed probe 
density, for example, a given microarray surface area 
of the defined subset genome-derived single exon 
microarray can yield a greater number of expression 
measurements. Alternatively, at a given probe density, 
the same number of expression measurements can be 
obtained from a smaller substrate surface area. 
Alternatively, at a fixed probe density and fixed 
surface area, probes can be provided redundantly, 
providing greater reliability in signal measurement for 
any given probe. Furthermore, with a higher percentage 
of probes known to be expressed in the assayed tissue, 
the dynamic range of the detection means can be 
adjusted to reveal finer levels discrimination among 
the levels of expression. 

In another aspect of the present invention, a 
genome-derived single-exon microarray is packaged 
together with an addressable set of individual probes, 
the set of individual probes including at least a 
subset of the probes on the microarray. In alternative 
embodiments, the ordered set of amplifiable probes is 
packaged separately from the genome-derived single exon 
microarray. 

In some embodiments, the microarray and/or 
ordered probe set are further packaged with recorded 
media that provide probe identification and addressing 
information, and that can additionally contain 
annotation information, such as gene expression data. 



Such recorded media can be packaged with the 
microarray, with the ordered probe set, or with both. 

If the microarray is constructed on a 
substrate that incorporates recordable media, such as 
is described in international patent application no. 
WO 98/12559, entitled "Spatially addressable 
combinatorial chemical arrays in CD-ROM format, " 
incorporated herein by reference in its entirety, then 
separate packaging of the genome-derived single exon 
microarray and the bioinf ormatic information is not 
required. 

Although the use of high density 
genome-derived microarrays on solid planar substrates 
is presently a preferred approach for the physical 
confirmation and characterization of the expression of 
sequences predicted to encode protein, other types of 
microarrays, as well as lower density macro arrays, can 
also be used. 

Experimental verification in process 400 of 
the function predicted from genomic sequence in process 
200 can be bioinf ormatic, rather than, or additional 
to, physical verification. 

Where the function desired to be identified 
is protein coding, the predicted exons can be compared 
bioinformatically to sequences known or suspected of 
being expressed. 

Thus, the sequences output from process 300 
(or process 200), can be used to query expression 
databases, such as EST databases, SNP ("single 
nucleotide polymorphism") databases, known cDNA and 
mRNA sequences, SAGE ("serial analysis of gene 
expression") databases, and more generalized sequence 
databases that allow query for expressed sequences. 
Such query can be done by any sequence query algorithm, 
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such as BLAST ("basic local alignment search tool") . 
The results of such query — including information on 
identical sequences and information on nonidentical 
sequences that have diffuse or focal regions of 
5 sequence homology to the query sequence — can then be 
passed directly to process 500, or used to inform 
analyses subsequently undertaken in process 200, 
process 300, or process 400. 

Experimental data, whether obtained by 

10 physical or bioinf ormatic assay in process 400, is 

passed to process 500 where it is usefully related to 
the sequence data itself, a process colloquially termed 
"annotation". Such annotation can be done using any 
technique that usefully relates the functional 

15 information to the sequence, as, for example, by 
incorporating the functional data into the record 
itself, by linking records in a hierarchical or 
relational database, by linking to external databases, 
or by a combination thereof. Such database techniques 

20 are well within the skill in the art. 

The annotated sequence data can be stored 
locally, uploaded to genomic sequence database 100, 
and/or displayed 800. 

The methods and apparatus of the present 

25 invention rapidly produce functional information from 
genomic sequence. We have, for example, used the 
methods and apparatus of the present invention to 
identify over 15,000 exons in human genomic sequence 
whose expression we have confirmed in at least one 

30 human tissue or cell type. Fully two-thirds of the 

exons belong to genes that were not then represented in 
existing public expression (EST, cDNA) databases. We 
have also used these single exon probes to identify 
alternative splice events in novel genes. 
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Coupled with the escalating pace at which 
sequence now accumulates, the ability rapidly to 
identify and confirm the function of regions of genomic 
DNA provided by the present invention produces a need 
for methods of displaying the information in meaningful 
ways. It is, therefore, another aspect of the present 
invention to provide means for displaying annotated 
sequence, and in particular for displaying sequence 
annotated according to the methods and apparatus of the 
present invention. Further, such display can be used 
as a preferred graphical user interface for electronic 
search, query, and analysis of such annotated sequence. 

FIG. 3 schematizes visual display 80 
presenting a single genomic sequence annotated 
according to the present invention. Because of its 
nominal resemblance to artistic works of Piet Mondrian, 
visual display 80 is alternatively described herein as 
a "Mondrian" . 

Each of the visual elements of display 80 is 
aligned with respect to the genomic sequence being 
annotated (the "annotated sequence") . Given the number 
of nucleotides typically represented in an annotated 
sequence, representation of individual nucleotides 
would rarely be readable in hard copy output of display 
80. Typically, therefore, the annotated sequence is 
schematized as rectangle 89, extending from the left 
border of display 80 to its right border. By 
convention herein, the left border of rectangle 89 
represents the first nucleotide of the sequence and the 
right border of rectangle 89 represents the last 
nucleotide of the sequence. 

As further discussed below, however, the 
Mondrian visual display of annotated sequence can serve 
as a convenient graphical user interface for 
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computerized representation, analysis, and query of 
information stored electronically. For such use, the 
individual nucleotides can conveniently be linked to 
the X axis coordinate of rectangle 89. This permits 
5 the annotated sequence at any point within rectangle 89 
readily to be viewed, either automatically — for 
example, by time-delayed appearance of a small overlaid 
window ("tool tip") upon movement of a cursor or other 
pointer over rectangle 89 — or through user 

10 intervention, as by clicking a mouse or other pointing 
device at a point in rectangle 89. 

Visual display 80 is generated after user 
specification of the genomic sequence to be displayed. 
Such specification can consist of or include an 

15 accession number for a single clone (e.g., a single BAC 
accessioned into GenBank) , wherein the starting and 
stopping nucleotides are thus absolutely identified, or 
alternatively can consist of or include an anchor or 
fulcrum point about which a chosen range of sequence is 

20 anchored, thus providing relative endpoints for the 
sequence to be displayed. For example, the user can 
anchor such a range about a given chromosomal map 
location, gene name, or even a sequence returned by 
query for similarity or identity to an input query 

25 sequence. When visual display 80 is used as a 
graphical user interface to computerized data, 
additional control over the first and last displayed 
nucleotide will typically be dynamically selectable, as 
by use of standard zooming and/or selection tools. 

30 Field 81 of visual display 80 is used to 

present the output from process 200, that is, to 
present the bioinf ormatic prediction of those sequences 
having the desired function within the genomic 
sequence. Functional sequences are typically indicated 
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by at least one rectangle 83 (83a, 83b, 83c), the left 
and right borders of which respectively indicate, by 
their X-axis coordinates, the starting and ending 
nucleotides of the region predicted to have function. 

Where a single bioinf ormatic method or 
approach identifies a plurality of regions having the 
desired function, a plurality of rectangles 83 is 
disposed horizontally in field 81. Where multiple 
methods and/or approaches are used to identify 
function, each such method and/or approach can be 
represented by its own series of horizontally disposed 
rectangles 83, each such horizontally disposed series 
of rectangles offset vertically from those representing 
the results of the other methods and approaches. 

Thus, rectangles 83a in FIG. 3 represent the 
functional predictions of a first method of a first 
approach for predicting function, rectangles 83b 
represent the functional predictions of a second method 
and/or second approach for predicting that function, 
and rectangles 83c represent the predictions of a third 
method and/or approach. 

Where the function desired to be identified 
is protein coding, field 81 is used to present the 
bioinformatic prediction of sequences encoding protein. 
For example, rectangles 83a can represent the results 
from GRAIL or GRAIL II, rectangles 83b can represent 
the results from GENE FINDER, and rectangles 83c can 
represent the results from DICTION. 

Optionally, and preferably, rectangles 83 
collectively representing predictions of a single 
method and/or approach are identically colored and/or 
textured, and are distinguishable from the color and/or 
texture used for a different method and/or approach. 



- 65 - 



Alternatively, or in addition, the color, 
hue, density, or texture of rectangles 83 can be used 
further to report a measure of the bioinf ormatic 
reliability of the prediction. For example, many gene 
prediction programs will report a measure of the 
reliability of prediction. Thus, increasing degrees of 
such reliability can be indicated, e.g., by increasing 
density of shading. Where display 80 is used as a 
graphical user interface, such measures of reliability, 
and indeed all other results output by the program, can 
additionally or alternatively be made accessible 
through linkage from individual rectangles 83, as by 
time-delayed window ("tool tip" window), or by pointer 
(e.g., mouse) -activated link. 

As above described, increased predictive 
reliability can be achieved by requiring consensus 
among methods and/or approaches to determining 
function. Thus, field 81 can include a horizontal 
series of rectangles 83 that indicate one or more 
degrees of consensus in predictions of function, 
including the combined length of the separately 
predicted exons that overlap in frame. 

Although FIG. 3 shows three series of 
horizontally disposed rectangles in field 81, 
display 80 can include as few as one such series of 
rectangles and as many as can discriminably be 
displayed, depending upon the number of methods and/or 
approaches used to predict a given function. For 
example, addition of a fourth gene prediction program, 
such as GENSCAN 

(http : //genes .mit . edu/GENSCANinf o . html ) , 

to the three gene prediction programs used in our first 
experiments (GRAIL, GENEFINDER, DICTION) would be 
accommodated by a fourth series of rectangles disposed 
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horizontally in field 81, but offset vertically from 
rectangles 81a, 81b, and 81c. 

Furthermore, field 81 can be used to show 
predictions of a plurality of different functions. 
5 However, the increased visual complexity occasioned by 
such display makes more useful the ability of the user 
to select a single function for display. When display 
80 is used as a graphical user interface for computer 
query and analysis, such function can usefully be 

10 indicated and user-selectable, as by a series of 
graphical buttons or tabs (not shown in FIG. 3) . 

Rectangle 89 is shown in FIG . 3 as including 
interposed rectangle 84. Rectangle 84 represents the 
portion of annotated sequence for which predicted 

15 functional information has been assayed physically, 

with the starting and ending nucleotides of the assayed 
material indicated by the X axis coordinates of the 
left and right borders of rectangle 84. Rectangle 85, 
with optional inclusive circles 86 (86a, 86b, and 86c) 

20 displays the results of such physical assay. 

Although a single rectangle 84 is shown in 
FIG. 3, physical assay is not limited to just one 
region of annotated genomic sequence. It is expected 
that an increasing percentage of regions predicted to 

25 have function by process 200 will be assayed 

physically, and that display 80 will accordingly, for 
any given genomic sequence, have an increasing number 
of rectangles 84 and 85, representing an increased 
density of sequence annotation. For example, for 

30 purposes of generating exon-specif ic probes for 

alternative splice detection, it is preferred that a 
plurality of exons, preferably all of the exons, that 
commonly belong to a single gene will be assayed 
experimentally for expression; accordingly, display 80 
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will have, for the genomic sequence encompassing such 
exons, a series of rectangles 84 and 85 for each of the 
assayed exons. 

Where the function desired to be identified 
5 is protein coding, rectangle 84 identifies the sequence 
of the probe used to measure expression. In 
embodiments of the present invention where expression 
is measured using genome-derived single exon 
microarrays, rectangle 84 identifies the sequence 

10 included within the probe immobilized on the solid 
support surface of the microarray. As noted supra, 
such probe will often include a small amount of 
additional, synthetic, material incorporated during 
amplification and designed to permit reamplif ication of 

15 the probe, which sequence is typically not shown in 
display 80. 

Rectangle 87 is used to present the results 
of bioinf ormatic assay of the genomic sequence. For 
example, where the function desired to be identified is 

20 protein coding, process 400 can include bioinf ormatic 
query of expression databases with the sequences 
predicted in process 200 to encode exons. And as above 
discussed, because bioinf ormatic assay presents fewer 
constraints than does physical assay, often the entire 

25 output of process 200 can be used for such assay, 
without further subsetting thereof by process 300. 
Therefore, rectangle 87 typically need not have 
separate indicators therein of regions submitted for 
bioinf ormatic assay; that is, rectangle 87 typically 

30 need not have regions therein analogous to 
rectangles 84 within rectangle 89. 

Rectangle 87 as shown in FIG. 3 includes 
smaller rectangles 880 and 88. Rectangles 880 indicate 
regions that returned a positive result in the 
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bioinformatic assay, with rectangles 88 representing 
regions that did not return such positive results. 
Where the function desired to be predicted and 
displayed is protein coding, rectangles 880 indicate 
5 regions of the predicted exons that identify sequence 
with significant similarity in expression databases, 
such as EST, SNP, SAGE databases, with rectangles 88 
indicating genes novel over those identified in 
existing expression data bases. 

10 Rectangles 880 can further indicate, through 

color, shading, texture, or the like, additional 
information obtained from bioinformatic assay. 

For example, where the function assayed and 
displayed is protein coding, the degree of shading of 

15 rectangles 880 can be used to represent the degree of 
sequence similarity found upon query of expression 
databases. The number of levels of discrimination can 
be as few as two (identity, and similarity, where 
similarity has a user-selectable lower threshold) . 

20 Alternatively, as many different levels of 

discrimination can be indicated as can visually be 
discriminated. 

Where display 80 is used as a graphical user 
interface, rectangles 880 can additionally provide 

25 links directly to the sequences identified by the query 
of expression databases, and/or statistical summaries 
thereof. As with each of the precedingly-discussed 
uses of display 80 as a graphical user interface, it 
should be understood that the information accessed via 

30 display 80 need not be resident on the computer 

presenting such display, which often will be serving as 
a client, with the linked information resident on one 
or more remotely located servers. 
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Rectangle 85 displays the results of physical 
assay of the sequence delimited by its left and right 
borders . 

Rectangle 85 can consist of a single 
5 rectangle, thus indicating a single assay, or 

alternatively, and increasingly typically, will consist 
of a series of rectangles (85a, 85b, 85c) indicating 
separate physical assays of the same sequence. 

Where the function assayed is gene 

10 expression, and where gene expression is assayed as 
herein described using simultaneous two-color 
fluorescent detection of hybridization to genome- 
derived single exon microarrays, individual rectangles 
85 can be colored to indicate the degree of expression 

15 relative to control. Conveniently, shades of green can 
be used to depict expression in the sample over control 
values, and shades of red used to depict expression 
less than control, corresponding to the spectra of the 
Cy3 and Cy5 dyes conventionally used for respective 

20 labeling thereof. Additional functional information 
can be provided in the form of circles 86 (86a, 86b, 
86c) , where the diameter of the circle can be used to 
indicate a parameter different from that set forth in 
rectangle 85. For example, where the annotated 

25 functions are the distribution of expression of the one 
or more predicted exons, rectangle 85 can report 
expression relative to control and circle 8 6 can be 
used to report signal intensity. As discussed infra, 
such relative expression (expression ratio) and 

30 absolute expression (signal intensity) can be expressed 
using normalized values. 

Where display 80 is used as a graphical user 
interface, rectangle 85 can be used as a link to 
further information about the assay. For example, 
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where the assay is one for gene expression, each 
rectangle 85 can be used to link to information about 
the source of the hybridized mRNA, the identity of the 
control, raw or processed data from the microarray 
5 scan, or the like. 

For purposes of illustration only, FIG. 4 
shows an embodiment of display 80 showing typical color 
conventions when hypothetical genomic sequence is 
annotated with exon-specif ic expression data. As would 

10 of course readily be understood, the color choice is 
arbitrary, and alternative colors can be used. 

In this typical presentation, BAC sequence 
("Chip seq.") 89 is presented in red, with the 
physically assayed region thereof (corresponding to 

15 rectangle 84 in FIG. 3) shown in white. Algorithmic 
gene predictions are shown in field 81, with 
predictions by GRAIL shown in green, predictions by 
GENEFINDER shown in blue, and predictions by DICTION 
shown in pink. Within rectangle 87, regions of 

20 sequence that, when used to query expression databases, 
return identical or similar sequences ("EST hit") are 
shown as white rectangles (corresponding to rectangles 
880 in FIG. 3) , gray indicates low homology, and black 
indicates unknowns (where black and gray would 

25 correspond to rectangles 88 in FIG. 3) . 

Although FIGS. 3 and 4 show a single stretch 
of sequence, uninterrupted from left to right, longer 
sequences are usefully represented by vertical stacking 
of such individual Mondrians, as shown in FIGS. 9 and 

30 10. 

Using our visual display tool, the Mondrian, 
we have found that consensus in the pattern of 
expression of individual exons is a powerful means for 
identifying exons that commonly belong to a single 
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gene. It is, therefore, another aspect of the present 
invention to provide methods, including methods based 
upon visual display, for associating exons that 
commonly belong to a single gene using, as the 
5 criterion for association, consensus in their patterns 
of expression in a plurality of tissues and/or cell 
types . 

As further discussed in Example 3, FIG. 9 
presents a Mondrian of BAC AC008172 (bases 25,000 to 

10 130,000 shown), containing the carbamyl phosphate 
synthetase gene (AF154830 . 1 ) , the sequence and 
structure of which has previously been reported. 
Purple background within the region shown as field 81 
in FIG. 3 indicates all 37 known exons for this gene. 

15 As can be seen, GRAIL II successfully 

identified 27 of the known exons (73%), GENE FINDER 
successfully identified 37 of the known exons (100%) , 
while DICTION identified 7 of the known exons (19%) . 

Seven of the predicted exons were selected 

20 for physical assay, of which 5 successfully amplified 
by PCR and were sequenced. These five exons were all 
found to be from the same gene, the carbamyl phosphate 
synthetase gene (AF154830 . 1 ) . 

The five exons were arrayed and gene 

25 expression measured across 10 tissues. As is readily 
seen by visual inspection of the resulting Mondrian 
(FIG. 5), the five single-exon probes report identical 
expression ratio patterns: each exon is expressed above 
control (i.e., in green) in the tissues represented by 

30 the fourth, seventh, and eighth rectangles 

(corresponding to rectangles 85 in FIG. 3) and is 
expressed at or below control in the remaining tissues. 

Of course, an exon that is removed or 
truncated by alternatively splicing in one of the 
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assayed tissues would produce a variant expression 
pattern. For purposes of associating exons as 
belonging commonly to a single gene, however, a 
consensus among assayed tissues would still identify 
5 the exon as presumptively belonging to the same gene. 

The methods of this aspect of the invention 
can, and typically will, be automated. For example, 
WO 99/58720, incorporated herein by reference in its 
entirety, describes algorithms for ordering the 

10 relatedness of a plurality of multidimensional 

expression data sets. The methods set forth therein 
can readily be adapted to ordering the relatedness of 
data sets, wherein each data set comprises expression 
ratios of an individual exon across a plurality of 

15 tissues and cell types, permitting exons with related, 
but not necessarily identical, patterns of expression 
to be classified as belonging to a common gene. 



The following examples are offered by way of 
20 illustration and not by way of limitation. 

EXAMPLE 1 

Preparation of Single Exon Microarrays 
from Exons Predicted in Human Genomic Sequence 

Bioinf ormatics Results 

25 All human BAC sequences in fewer than 10 

pieces that had been accessioned in a five month period 
immediately preceding this study were downloaded from 
GenBank. This corresponds to -2200 clones, totaling 
-350 MB of sequence, or approximately 10% of the human 

3 0 genome. 
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After masking repetitive elements using the 
program CROSS_MATCH, the sequence was analyzed for open 
reading frames using three separate gene finding 
programs. The three programs predict genes using 
5 independent algorithmic methods developed on 

independent training sets: GRAIL uses a neural network, 
GENE FINDER uses a hidden Markoff model, and DICTION, a 
program proprietary to Genetics Institute, operates 
according to a different heuristic. The results of all 

10 three programs were used to create a prediction matrix 
across the segment of genomic DNA. 

The three gene finding programs yielded a 
range of results. GRAIL identified the greatest 
percentage of genomic sequence as putative coding 

15 region, 2% of the data analyzed. GENEFINDER was 
second, calling 1%, and DICTION yielded the least 
putative coding region, with 0.8% of genomic sequence 
called as coding region. 

The consensus data were as follows. GRAIL 

20 and GENEFINDER agreed on 0.7% of genomic sequence, 

GRAIL and DICTION agreed on 0.5% of genomic sequence, 
and the three programs together agreed on 0.25% of the 
data analyzed. That is, 0.25% of the genomic sequence 
was identified by all three of the programs as 

25 containing putative coding region. 

Exons predicted by any two of the three 
programs ("consensus exons") were assorted into "gene 
bins" using two criteria: (1) any 7 consecutive exons 
within a 25 kb window were placed together in a bin as 

30 likely contributing to a single gene, and (2) all exons 
within a 25 kb window were placed together in a bin as 
likely contributing to a single gene if fewer than 7 
exons were found within the 25 kb window. 
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PCR 

The largest exon from each gene bin that did 
not span repetitive sequence was then chosen for 
amplification/ as were all consensus exons longer than 
5 500 bp. This method approximated one exon per gene; 

however^- a number of genes were found to be represented 
by multiple elements. 

Previously, we had determined that DNA 
fragments fewer than 250 bp in length do not bind well 

10 to the amino-modif ied glass surface of the slides used 
as support substrate for construction of microarrays; 
therefore, amplicons were designed in the present 
experiments to approximate 500 bp in length. 

Accordingly, after selecting the largest exon 

15 per gene bin, a 500 bp fragment of sequence centered on 
the exon was passed to the primer picking software, 
PRIMER3 (available online for use at 
http://www-genome.wi.mit.edu/cgi-bin/primer/ ) . A 
first additional sequence was commonly added to each 

20 exon-unique 5 1 primer, and a second, different, 

additional sequence was commonly added to each exon- 
unique 3 T primer, to permit subsequent reamplif ication 
of the amplicon using a single set of "universal" 5 T 
and 3 1 primers, thus immortalizing the amplicon. The 

25 addition of universal priming sequences also 

facilitates sequence verification, and can be used to 
add a cloning site should some exons be found to 
warrant further study. 

The exons were then PCR amplified from 

30 genomic DNA, verified on agarose gels, and sequenced 

using the universal primers to validate the identity of 
the amplicon to be spotted in the microarray. 

Primers were supplied by Operon Technologies 
(Alameda, CA) . PCR amplification was performed by 



standard techniques using human genomic DNA (Clontech, 
Palo Alto, CA) as template. Each PCR product was 
verified by SYBR® green (Molecular Probes, Inc., 
Eugene, OR) staining of agarose gels, with subsequent 
imaging by Fluorimager (Molecular Dynamics, Inc., 
Sunnyvale, CA) . PCR amplification was classified as 
successful if a single band appeared. 

The success rate for amplifying exons of 
interest directly from genomic DNA using PCR was 
approximately 75%. FIG. 5 graphs the distribution of 
predicted exon length and distribution of amplified PCR 
products, with exon length shown by dashed line and PCR 
product length shown by solid line. Although the range 
of exon sizes is readily seen to extend to beyond 
900 bp, the mean predicted exon size was only 229 bp, 
with a median size of 150 bp (n-9498) . With an average 
amplicon size of 475 ± 25 bp, approximately 50% of the 
average PCR amplification product contained predicted 
coding region, with the remaining 50% of the amplicon 
containing either intron, intergenic sequence, or both. 

Using a strategy predicated on amplifying 
about 500 bp, it was found that long exons had a higher 
PCR failure rate. To address this, the bioinf ormatics 
process was adjusted to amplify 1000, 1500 or 2000 bp 
fragments from exons larger than 500 bp. This improved 
the rate of successful amplification of exons exceeding 
500 bp, constituting about 9.2% of the exons predicted 
by the gene finding algorithms. 

Approximately 75% of the probes disposed on 
the array (90% of those that successfully PCR 
amplified) were sequence-verified by sequencing in both 
the forward and reverse direction using MegaBACE 
sequencer (Molecular Dynamics, Inc., Sunnyvale, CA) , 
universal primers, and standard protocols. 



Some genomic clones (BACs) yielded very poor 
PCR and sequencing results. The reasons for this are 
unclear, but may be related to the quality of early 
draft sequence or the inclusion of vector and host 
contamination in some submitted sequence data. 

Although the intronic and intergenic material 
flanking coding regions could theoretically interfere 
with hybridization during microarray experiments, 
subsequent empirical results demonstrated that 
differential expression ratios were not significantly 
affected by the presence of noncoding sequence. The 
variation in exon size was similarly found not to 
affect differential expression ratios significantly; 
however, variation in exon size was observed to affect 
the absolute signal intensity (data not shown) . 

The 350 MB of genomic DNA was, by the above- 
described process, reduced to 9750 discrete probes, 
which were spotted in duplicate onto glass slides using 
commercially available instrumentation (MicroArray 
Genii Spotter and/or MicroArray Genlll Spotter, 
Molecular Dynamics, Inc., Sunnyvale, CA) . Each slide 
additionally included either 16 or 32 E. coli genes, 
the average hybridization signal of which was used as a 
measure of background biological noise. 

Each of the probe sequences was BLASTed 
against the human EST data set, the NR data set, and 
SwissProt GenBank (May 7, 1999 release 2.0.9). 

One third of the probe sequences (as 
amplified) produced an exact match (BLAST Expect ("E") 
values less than 1 e" 100 ) to either an EST (20% of 
sequences) or a known mRNA (13% of sequences) . A 
further 22% of the probe sequences showed some homology 
to a known EST or mRNA (BLAST E values from 1 e~ 5 to 
1 e~") . The remaining 45% of the probe sequences 



showed no significant sequence homology to any 
expressed, or potentially expressed, sequences present 
in public databases. 

All of the probe sequences (as amplified) 
were then analyzed for protein similarities with the 
SwissProt database using BLASTX, Gish et al . f Nature 
Genet. 3:266 (1993). The predicted functional 
breakdowns of the 2/3 of probes identical or homologous 
to known sequences are presented in Table 1. 



Table 1 



Function of Predicted Exons 
As Deduced From Comparative Sequence Analysis 

Total V6 chip V7 chip Function Predicted from 

Comparative Sequence 
Analysis 


211 


96 


115 


Receptor 


120 


43 


77 


Zinc Finger 


30 


11 


19 


Homeobox 


25 


9 


16 


Transcription Factor 


17 


11 


7 


Transcription 


118 


57 


61 


Structural 


95 


39 


56 


Kinase 


36 


18 


18 


Phosphatase 


83 


31 


52 


Ribosomal 


45 


19 


26 


Transport 


21 


7 


14 


Growth Factor 


17 


12 


5 


Cytochrome 


50 


33 


17 


Channel 



As can be seen, the two most common types of 
genes were transcription factors and receptors, making 
up 2.2% and 1.8% of the arrayed elements, respectively. 
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EXAMPLE 2 
Gene Expression Measurements From 
Genome-Derived Single Exon Microarrays 



The two genome-derived single exon 
5 microarrays prepared according to Example 1 were 
hybridized in a series of simultaneous two-color 
fluorescence experiments to (1) Cy3-labeled cDNA 
synthesized from message drawn individually from each 
of brain, heart, liver, fetal liver, placenta, lung, 

10 bone marrow, HeLa, BT 474, or HBL 100 cells, and (2) 

Cy5-labeled cDNA prepared from message pooled from all 
ten tissues and cell types, as a control in each of the 
measurements. Hybridization and scanning were carried 
out using standard protocols and Molecular Dynamics 

15 equipment. 

Briefly, mRNA samples were bought from 
commercial sources (Clontech, Palo Alto, CA and 
Amersham Pharmacia Biotech (APB) ) . Cy3-dCTP and 
Cy5-dCTP (both from APB) were incorporated during 

20 separate reverse transcriptions of 1 pg of polyA + mRNA 
performed using 1 pg oligo (dT) 12-18 primer and 2 jig 
random 9mer primers as follows. After heating to 70°C, 
the RNA: primer mixture was snap cooled on ice. After 
snap cooling on ice, added to the RNA to the stated 

25 final concentration was: IX Superscript II buffer, 0.01 
M DTT, 100]jM dATP, 100 pM dGTP, 100 pM dTTP, 50 pM 
dCTP, 50 pM Cy3-dCTP or Cy5-dCTP 50 pM, and 200 U 
Superscript II enzyme. The reaction was incubated for 
2 hours at 42°C. After 2 hours, the first strand cDNA 

30 was isolated by adding 1 U Ribonuclease H, and 

incubating for 30 minutes at 37°C. The reaction was 
then purified using a Qiagen PCR cleanup column, 
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increasing the number of ethanol washes to 5. Probe 
was eluted using 10 mM Tris pH 8.5. 

Using a spectrophotometer, probes were 
measured for dye incorporation. Volumes of both Cy3 
5 and Cy5 cDNA corresponding to 50 pmoles of each dye 
were then dried in a Speedvac, resuspended in 30 \il 
hybridization solution containing 50% formamide, 
5X SSC, 0.2 ug/ul poly(dA), 0.2 jag/pl human c 0 tl DNA, 
and 0.5 % SDS. 

10 Hybridizations were carried out under a 

coverslip, with the array placed in a humid oven at 
42°C overnight. Before scanning, slides were washed in 
IX SSC, 0.2% SDS at 55°C for 5 minutes, followed by 
0.1X SSC, 0.2% SDS, at 55°C for 20 minutes. Slides 

15 were briefly dipped in water and dried thoroughly under 
a gentle stream of nitrogen. 

Slides were scanned using a Molecular 
Dynamics Gen3 scanner, as described. Schena (ed.), 
Microarray Biochip: Tools and Technology / Eaton 

20 Publishing Company/BioTechniques Books Division (2000) 
(ISBN; 1881299376) . 

Although the use of pooled cDNA as a 
reference permitted the survey of a large number of 
tissues, it attenuates the measurement of relative gene 

25 expression, since every highly expressed gene in the 
tissue/cell type-specific fluorescence channel will be 
present to a level of at least 10% in the control 
channel. Because of this fact, both signal and 
expression ratios (the latter hereinafter, "expression" 

30 or "relative expression") for each probe were 

normalized using the average ratio or average signal, 
respectively, as measured across the whole slide. 

Data were accepted for further analysis only 
when signal was at least three times greater than 
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biological noise, the latter defined by the average 
signal produced by the E. coli control genes. 

The relative expression signal for these 
probes was then plotted as a function of tissue or cell 
5 type, and is presented in FIG . 6. 

FIG. 6 shows the distribution of expression 
across a panel of ten tissues. The graph shows the 
number of sequence-verified products that were either 
not expressed ("0"), expressed in one or more but not 
10 all tested tissues ("1" - "9"), and expressed in all 
tissues tested ("10") . 

Of 9999 arrayed elements on the two 
microarrays (including positive and negative controls 
and "failed" products), 2353 (51%) were expressed in at 
15 least one tissue or cell type. Of the gene elements 
showing significant signal — where expression was 
scored as "significant" if the normalized Cy3 signal 
was greater than 1, representing signal 5-fold over 
biological noise (0.2) — 39% (991) were expressed in 
20 all 10 tissues. The next most common class (15%) 

consisted of gene elements expressed in only a single 
tissue . 

The genes expressed in a single tissue were 
further analyzed, and the results of the analyses are 

25 compiled in FIG. 7. 

FIG. 7A is a matrix presenting the expression 
of all verified sequences that showed signal intensity 
greater than 3 in at least one tissue. Each clone is 
represented by a column in the matrix. Each of the 10 

30 tissues assayed is represented by a separate row in the 
matrix, and relative expression (expression ratio) of a 
clone in that tissue is indicated at the respective 
node by intensity of green shading, with the intensity 
legend shown in panel B. The top row of the matrix 
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("EST Hit") contains "bioinf ormatic" rather than 
"physical" expression data — that is, presents the 
results returned by query of EST, NR and SwissProt 
databases using the probe sequence. The legend for 
5 "bioinf ormatic expression" {i.e., degree of homology 
returned) is presented in panel C. Briefly, white is 
known, black is novel, with gray depicting nonidentical 
with significant homology (white: E values < 1 e~ 100 ; 
gray: E values from le" 5 (1 x 10" 5 ) to le~" (1 x 10""); 

10 black: E values > le" 5 (1 x 10" 5 ) . 

As FIG . 7 readily shows, heart and brain were 
demonstrated to have the greatest numbers of genes that 
were shown to be uniquely expressed in the respective 
tissue. In brain, 200 uniquely expressed genes were 

15 identified; in heart, 150. The remaining tissues gave 
the following figures for uniquely expressed genes: 
liver, 100; lung, 70; fetal liver, 150; bone marrow, 
75; placenta, 100; HeLa, 50; HBL, 100; and BT474, 50. 

It was further observed that there were many 

20 more "novel" genes among those that were up-regulated 
in only one tissue, as compared with those that were 
down-regulated in only one tissue. In fact, it was 
found that exons whose expression was measurable in 
only a single of the tested tissues were represented in 

25 sequencing databases at a rate of only 11%, whereas 36% 
of the exons whose expression was measurable in 9 of 
the tissues were present in public databases. As for 
those exons expressed in all ten tissues, fully 45% 
were present in existing expressed sequence databases. 

30 These results are not unexpected, since genes expressed 
in a greater number of tissues have a higher likelihood 
of being, and thus of having been, discovered by EST 
approaches . 
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Comparison of Signal from Known and Unknown Genes 

The normalized signal of the genes found to 
have high homology to genes present in the GenBank 
human EST database were compared to the normalized 
5 signal of those genes not found in the GenBank human 
EST database. The data are shown in FIG. 8. 

FIG. 8 shows in dashed line the normalized 
Cy3 signal intensity for all sequence-verified products 
with a BLAST Expect ("E") value of greater than le~ 30 

10 (1 x 10" 30 ) (designated "unknown") upon query of 

existing EST, NR and SwissProt databases, and shows in 
blue the normalized Cy3 signal intensity for all 
sequence-verified products with a BLAST Expect value of 
less than le" 30 ("known"). Note that biological 

15 background noise has an averaged normalized Cy3 signal 
intensity of 0.2. 

As expected, the most highly expressed of the 
exons were "known" genes. This is not surprising, 
since very high signal intensity correlates with very 

20 commonly-expressed genes, which have a higher 
likelihood of being found by EST sequence. 

However, a significant point is that a large 
number of even the high expressers were "unknown" . 
Since the genomic approach used to identify genes and 

25 to confirm their expression does not bias exons toward 
either the 3 T or 5 T end of a gene, many of these high 
expression genes will not have been detected in an end- 
sequenced cDNA library. 

The significant point is that presence of the 

30 gene in an EST database is not a prerequisite for 

incorporation into a genome-derived microarray, and 
further, that arraying such "unknown" exons can help to 
assign function to as-yet undiscovered genes. 
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Verification of Gene Expression 

To ascertain the validity of the approach 
described above to identify genes from raw genomic 
sequence, expression of two of the probes was assayed 
5 using reverse transcriptase polymerase chain reaction 
(RT PCR) and northern blot analysis. 

Two microarray probes were selected on the 
basis of exon size, prior sequencing success, and 
tissue-specific gene expression patterns as measured by 
10 the microarray experiments. The primers originally 
used to amplify the two respective exons from genomic 
DNA were used in RT PCR against a panel of tissue- 
specific cDNAs (Rapid-Scan gene expression panel 24 
human cDNAs) (OriGene Technologies, Inc., Rockville, 
15 MD) . 

Sequence AL079300_1 was shown by microarray 
hybridization to be present in cardiac tissue, and 
sequence AL031734_1 was shown by microarray experiment 
to be present in placental tissue (data not shown) . 

20 RT-PCR on these two sequences confirmed the tissue- 
specific gene expression as measured by microarrays, as 
ascertained by the presence of a correctly sized PCR 
product from the respective tissue type cDNAs . 

Clearly, all microarray results cannot, and 

25 indeed should not, be confirmed by independent assay 
methods, or the high throughput, highly parallel 
advantages of microarray hybridization assays will be 
lost. However, in addition to the two RT-PCR results 
presented above, the observation that 1/3 of the 

30 arrayed genes exist in expression databases provides 

powerful confirmation of the power of our methodology — 
which combines bioinf ormatic prediction with expression 
confirmation using genome-derived single exon 
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microarrays — to identify novel genes from raw genomic 
data . 

To verify that the approach further provides 
correct characterization of the expression patterns of 
5 the identified genes, a detailed analysis was performed 
of the microarrayed sequences that showed high signal 
in brain. 

For this latter analysis, sequences that 
showed high (normalized) signal in brain, but which 

10 showed very low (normalized) signal (less than 0.5, 
determined to be biological noise) in all other 
tissues, were further studied. There were 82 sequences 
that fit these criteria, approximately 2% of the 
arrayed elements. The 10 sequences showing the highest 

15 signal in brain in microarray hybridizations are 

detailed in Table 2, along with assigned function, if 
known or reasonably predicted. 



Table 2 







Function of the 


Most Highly 


20 


Expressed Genes Expressed Only 


in Brain 




Microarray 


Normal 


Express 


Homology 


Gene Function 




Sequence 


ized 


ion 


to EST 


as described 




Name 


Signal 


Ratio 


present 


by GenBank 










in 












GenBank 






AP000217-1 


5.2 


+ 7.7 


High 


S-100 protein, 












b-chain, Ca 2+ 












binding 












protein 












expressed in 












central 












nervous system 


25 


AP000047-1 


2.3 




High 


Unknown 












Function 




AC006548-9 


1.7 




High 


Similar to 












mouse membrane 












glyco -protein 












M6, expressed 
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Function of the Most Highly 
Expressed Genes Expressed Only in Brain 

Microarray Normal Express Homology Gene Function 
Sequence ized ion to EST as described 
Name Signal Ratio present by GenBank 

in 

GenBank 










in central 
nervous system 


AC007245-5 


1.5 




High 


Similar to 
amphiphysin, a 
synapt i c 
vesicle- 
associated 
protein . Ref 
21 


L44140-4 


1.2 


+ 2.0 


High 


Endothelial 
act in-binding 
protein found 
in nonmuscle 
f ilamin 


AC004689-9 


1.2 


+ 3.5 


High 


Protein 

Phosphatase 

PP2A, 

neuronal/ 

downregulates 

activated 

protein 

kinases 


AL031657-1 


1 . 2 


+ 3.0 


High 


Unknown 
function/ 

r'^ri +■ 2 i "nc -f~ Vi o 
^OilLcLXXlo Ullti 

anhyrin motif, 
a common 
protein 
sequence motif 


AC009266-2 


1.1 


+ 3.7 


Low 


Low homology 
to the 

Synaptotagmin 
I protein in 
rat/present at 
low levels 
throughout rat 
brain 


AP000086-1 


1 . 0 


+ 2.7 


Low 


Unknown, very 
poor homology 
to collagen 


AC004689-3 


1.0 




High 


Protein 
Phosphatase 
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Function of the 


Most Highly 


Expressed Genes Expressed Only 


in Brain 


Microarray 


Normal 


Express 


Homology 




Sequence 


ized 


ion 


to EST 


as described 


Name 


Signal 


Ratio 


present 


by GenBank 








in 










GenBank 












PP2A, 










neuronal/ 










downregulates 










activated 










protein 










kinases 



Of the ten sequences studied by these latter 

35 confirmatory approaches, eight were previously known. 
Of these eight, six had previously been reported to be 
important in the central nervous system or brain. The 
exon giving the highest signal (AP00217-1) was found to 
be the gene encoding an S100B Ca 2+ binding protein, 

40 reported in the literature to be highly and uniquely 
expressed in the central nervous system. Heizmann, 
Neurochem. Res. 9:1097 (1997). 

A number of the brain-specific probe 
sequences (including AC006548-9, AC009266-2) did not 

4 5 have homology to any known human cDNAs in GenBank but 
did show homology to rat and mouse cDNAs . Sequences 
AC004689-9 and AC004689-3 were both found to be 
phosphatases present in neurons (Millward et al . , 
Trends Biochem. Sci. 24 ( 5) : 186-191 (1999)). Two 

50 microarray sequences, AP000047-1 and AP000086-1 have 
unknown function, with AP000086-1 being absent from 
GenBank. Functionality can now be narrowed down to a 
role in the central nervous system for both of these 
genes, showing the power of designing microarrays in 

55 this fashion. 
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Next, the function of the chip sequences with 
the highest (normalized) signal intensity in brain, 
regardless of expression in other tissues, was 
assessed. In this latter analysis, we found expression 
5 of many more common genes, since the sequences were not 
limited to those expressed only in brain. For example, 
looking at the 20 highest signal intensity spots in 
brain, 4 were similar to tubulin (AC00807905; AF146191- 
2; AC007664-4; AF14191-2), 2 were similar to actin 

10 (AL035701-2; AL034402-1), and 6 were found to be 

homologous to glycer aldehyde- 3 -phosphate dehydrogenase 
(GAPDH) (AL035604-1; Z86090-1; AC006064-L, AC006064-K; 
AC035604-3; AC006064-L) . These genes are often used as 
controls or housekeeping genes in microarray 

15 experiments of all types. 

Other interesting genes highly expressed in 
brain were a ferritin heavy chain protein, which is 
reported in the literature to be found in brain and 
liver (Joshi et al . , J. Neurol. Sci. 134 (Suppl) : 52-56 

20 (1995) ) , a result confirmed with the array. Other 

highly expressed chip sequences included a translation 
elongation factor la (AC007564-4 ) , a DEAD-box homolog 
(AL023804-4) , and a Y-chromosome RNA-binding motif 
(Chai et al . , Genomics 49(2):283-89 ( 1998 ) ) (AC007320- 

25 3) . A low homology analog (AP00123-1/2 ) to a gene, 
DSCR1, thought to be involved in trisomy 21 (Down's 
syndrome) , showed high expression in both brain and 
heart, in agreement with the literature (Fuentes et 
al., Mol. Genet. 4 ( 10) : 1935-44 (1995)). 

30 As a further validation of the approach, we 

selected the BAC AC006064 to be included on the array. 
This BAC was known to contain the GAPDH gene, and thus 
could be used as a control for the exon selection 
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process. The gene finding and exon selection 
algorithms resulted in choosing 25 exons from BAC 
AC006064 for spotting onto the array, of which four 
were drawn from the GAPDH gene* Table 3 shows the 
5 comparison of the average expression ratio for the 4 
exons from BAC0 0 60 64 compared with the average 
expression ratio for 5 different dilutions of a 
commercially available GAPDH cDNA (Clontech) . 



Table 3 



Comparison 


of Expression Ratio, 


for each tissue, 


of 




GAPDH 








AC006064 (n = 


4) Control ( n = 


5) 


Bone Marrow 


-1.81 ± 0.11 


-1.85 ± 0.08 




Brain 


-1.41 ± 0.11 


-1.17 ± 0.05 




BT474 


1.85 ± 0.09 


1.66 ± 0.12 




Fetal Liver 


-1.62 ± 0.07 


-1.41 ± 0.05 




HBL100 


1.32 ± 0.05 


2.64 ± 0.12 




Heart 


1.16 ± 0.09 


1.56 ± 0.10 




HeLa 


1.11 ±0.06 


1.30 ± 0.15 




Liver 


-1.62 ± 0.22 


-2.07 ± 




Lung 


-4.95 ± 0.93 


-3.75 ± 0.21 




Placenta 


-3.56 ± 0.25 


-3.52 ± 0.43 





Each tissue shows excellent agreement between 
the experimentally chosen exons and the control, again 
demonstrating the validity of the present exon mining 
25 approach. In addition, the data also show the 

variability of expression of GAPDH within tissues, 
calling into question its classification as a 
housekeeping gene and utility as a housekeeping control 
in microarray experiments. 
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EXAMPLE 3 
Representation of Sequence and 
Expression Data as a "Mondrian" 

For each genomic clone processed for 
5 microarray as above-described, a plethora of 

information was accumulated, including full clone 
sequence, probe sequence within the clone, results of 
each of the three gene finding programs, EST 
information associated with the probe sequences, and 
10 microarray signal and expression for multiple tissues, 
challenging our ability to display the information. 

Accordingly, we devised a new tool for visual 
display of the sequence with its attendant annotation 
which, in deference to its visual similarity to the 
15 paintings of Piet Mondrian, is hereinafter termed a 
"Mondrian". FIGS. 3 and 4 present the key to the 
information presented on a Mondrian. 

FIG. 9 presents a Mondrian of BAC AC008172 
(bases 25,000 to 130,000 shown), containing the 
20 carbamyl phosphate synthetase gene (AF154830 . 1 ) . 

Purple background within the region shown as field 81 
in FIG. 3 indicates all 37 known exons for this gene. 

As can be seen, GRAIL II successfully 
identified 27 of the known exons (73%), GENEFINDER 
25 successfully identified 37 of the known exons (100%), 
while DICTION identified 7 of the known exons (19%) . 

Seven of the predicted exons were selected 
for physical assay, of which 5 successfully amplified 
by PCR and were sequenced. These five exons were all 
30 found to be from the same gene, the carbamyl phosphate 
synthetase gene (AF154830 . 1 ) . 

The five exons were arrayed, and gene 
expression measured across 10 tissues. As is readily 
seen in the Mondrian, the five chip sequences on the 
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array show identical expression patterns, elegantly 
demonstrating the reproducibility of the system. 

FIG. 10 is a Mondrian of BAC AL049839. We 
selected 12 exons from this BAC , of which 10 
5 successfully sequenced, which were found to form 

between 5 and 6 genes. Interestingly, 4 of the genes 
on this BAC are protease inhibitors. Again, these data 
elegantly show that exons selected from the same gene 
show the same expression patterns, depicted below the 

10 red line. From this figure, it is clear that our 

ability to find known genes is very good. A novel gene 
is also found from 8 6.6 kb to 88.6 kb, upon which all 
the exon finding programs agree. We are confident we 
have two exons from a single gene since they show the 

15 same expression patterns and the exons are proximal to 
each other. Backgrounds in the following colors 
indicate a known gene (top to bottom) : 
red = kallistatin protease inhibitor (P29622) ; 
purple = plasma serine protease inhibitor (P05154); 

20 turquoise = al anti-chymotrypsin (P01011); mauve = 40S 
ribosomal protein (P08865) . Note that chip sequence 8 
and 12 did not sequence verify. 

EXAMPLE 4 
Sequences of Genes Identified From 
25 Genomic Sequence By Gene Prediction and Single Exon 

Microarray Analysis 



The sequences of three exons identified from 
human genomic sequence in experiments as set forth in 
Examples 1-3 are presented here, with each exon 
30 represented by its predicted coding sequence, and 

thereafter by the sequence of the amplicon as used on 
the genome-derived single exon microarray to assess its 
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expression. The three sequences were chosen, 
respectively, to represent each of three classes of 
genes obtainable by this method: (1) those that have 
already been identified and accessioned into expression 
5 databases such as EST, SNP, SwissProt databases; 
(2) those that are not identically represented in 
expression databases, but that have sequence showing 
significant homology to genes already present in such 
expression databases; and (3) those that are neither 
10 identically present nor have significant sequence 
homology to genes present in expression databases. 

The first, designated AC007683_4_chip. seq. 1, 
was found to be identical to a sequence in an existing 
expression database. 

15 ACQ07 683__4__ch ip.seq .l predicted exon : 

TTTTTTTTTTTGCAAGCAGATAAAGGCTTATTTTACTTTAATGGCTGATCTATGT 
AATCACGGAGGCCAGTATGTACACACAAAGGGGCAGCTTTTATTTCTTGGTCTCT 
TCCTCCTTGGACAAAGTCTTGATGATCTCCTCCTTCTTGGCCTGGAGGTGCTCTT 
CATAGCTCTTGTGTGCTTCCTTGGTCTTAGATCTGCGGGCCTCAGCCTGATCAGC 
2 0 CAGGAGCTTCTTGCGGGCCTTGTCTGCCTTCAGCTTGTGGATGTGTTCCATGAGA 
ATCTGCTTGTTTTTTAACACATTCCTCTTCACCTTCAGGTACAGGCTGTGATACA 
TGCGGCGATCAATCTTCTTA [SEQ ID NO:l] 

AC0 07 68 3_4_chip. seq. 1 amplicon : 

C AG T C C AC AT G G G T AC AA.G C C C T G AAAC C T C AAAT G T AC A T C AGAAT T AC C T G T G 
25 GAGTTGTTTTTTTTTTTTTTTTTTTTTTTTTGCAAGCAGATAAAGGCTTATTTTA 
C T T T AAT G G C T GAT C TAT G T AAT C AC G GAG GC C AG TAT G T AC AC AC AAAG G G G C A 
GCTTTTATTTCTTGGTCTCTTCCTCCTTGGACAAAGTCTTGATGATCTCCTCCTT 
CTTGGCCTGGAGGTGCTCTTCATAGCTCTTGTGTGCTTCCTTGGTCTTAGATCTG 
CGGGCCTCAGCCTGATCAGCCAGGAGCTTCTTGCGGGCCTTGTCTGCCTTCAGCT 
30 TGTGGATGTGTTCCATGAGAATCTGCTTGTTTTTTAACACATTCCTCTTCACCTT 
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CAGGTACAGGCTGTGATACATGCGGCGATCAATCTTCTTAGATTCACGGTATCTT 
CTGAGCAGCCGGTGCAGAATCCTCATTCTCCTCATCCACGTGACCTTCTCTGGCA 
TTCGG [SEQ ID NO : 2 ] . 

The second, designated AC007 682_2_chip . seq. 2, 
5 was not found identically in an expression database, 
but was found to have homology to one or more sequences 
in such databases. 

AC007 682_2_chip.seg.2 predicted exon : 

TAT GGTAT TTTCT T AT AGC AAC AAAAAAT AAAGAT GGGGTG GAGAAA T AT A T T TA 
10 T AGAAAGTAT T T T T T TAAGT [SEQ ID NO: 3] 

AC007 682_2_chi p. seer. 2 amplicon : 

AGTATGGAGCCCCCTTCATGGGACAGGTGGCTTTAAGAAGAGGAAGAGAGACCTG 
AGCTGGCAGGGACTCTCTTACCCTCTCACCATGTGATGCCCTCCACATGTTATGA 
TGCAGCAAGAAGGCCCTCACTGGTTGCTAGTGCCATGCTCTTCGACTTCCCAGCC 

1 5 T GCAGAAC T AT AAGAAAT AAACT T AT T T T C T T TAT AAC T T ACAC AT T TAT GGTAT 
TTTCT TAT AG CAACAAAAAAT AAAGAT GGGG T GGAGAAAT AT AT T T AT AGAAAG T 
AT T T T T T TAAG TAAAT GAGAAAT T AG AC AT AAT GTTTT T AAC T C T AGAGAAAT T G 
AAAAC AGAGCACAGCACAT C GGAT AAAT T C AAT AAC TAT C T T AAGAAT CAG CAAA 
ACAACAT GCAGATGGC TGAT TGGCAATAGT T T CAGT AGGCAGAT T T T GAT TAAAA 

2 0 T AAAGAAAAAC T T T T T AAT AAT TAAAC C T C T C C T T AAAAC AT TAT GAC T T TAT GA 
GGTAA [SEQ ID NO: 4] 

The third exon, designated 
AC007552_4_chip . seq. 2, was neither identically present 
nor significantly related in sequence to any entry in a 
25 public expression database. 
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AC0Q7552_4_chiio. seg.2 Predicted Exon 

TCTTCATTATTAATCACTCTTAZVACCTCTTCTTCAATCTTCTCCTCATGTTTAAT 
TTCTCCCTTATCTTATCTTCATAACTCAGTGCCATTCTCCCTTCATAACAACAGA 
AG C T GAC AT T GGAGG [SEQ ID NO: 5] 

5 ACQ07552_4__chip . sea. 2 amplicon : 

TCATCC TAAT T TATATAAAGCACAC TACAAT CTTAAT T T AAC AAT C CAT TCCAAA 
T T C CAAT AAT C T C C AG T G T T G AGAT AT T T T T T C CAT AC AG C C T AAAG T G C AC AT A 
TTTAGACATTTCTCCACCCATCTCCTTTGCACACGAAAAGTTGGTAAACGACCTC 
ATTATACTAGTAGCCTTTCATATTCTTCAT TATTAATCACTCTTAAACCTCTTCT 

10 TCAATCTTCTCCTCATGTTTAATTTCTCCCTTATCTTATCTTCATAACTCAGTGC 
CAT TCTCCCTT CAT AAC AAC AG AAG C T GAC AT T G GAG GAG TAT C AG C CAAT G T G T 
ACCGCTCTTTCCCTACTGTGGTCCACTGTCACCCCTAACTATTTTATGAATAGGA 
T T C C TAT T T C T AGAGAAGAAAAC GC AGAC T T G GAGAGG T T GAG T AAG T T GC C TAG 
GAATGTGAAGCTGGGGTGTAGCAGAAGGGGGTCGACGTCAGGTCTGGATACCTCA 

15 CCGTG [SEQ ID NO: 6] 



EXAMPLE 5 
Genome-Derived Single Exon Probes 
Useful For Measuring Human Gene Expression 

The protocols set forth in Examples 1 and 2, 
20 supra, were applied with some modification to 

additional human genomic sequence as it became newly 
available in GenBank. From the collective efforts of 
these and the experiments reported in Example 2, we 
generated over 15,000 unique human genome-derived 
25 single exon probes that could be shown to be expressed 
at significant levels in one or more of ten tested 
tissues . 
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Modifications to the protocols for 
bioinf ormatic prediction of exons set forth in Examples 
1 and 2 were as follows. 

First, we added a fourth gene prediction 
5 program, GENS CAN, to the three originally used, 
DICTION, GENEFINDER, and GRAIL . 

Second, we increased the resolution of our 
exon predictions, as follows. 

In the experiments reported In Examples 1 

10 and 2, we applied a 25 bp window in scanning genomic 
sequence: exons were called when any two of the three 
gene prediction programs identified an exon anywhere 
within the window. In the more recent experiments, we 
looked for consensus on a nucleotide by nucleotide 

15 basis: when any two or more of the four programs 

identified the nucleotide as falling within an exon, 
the nucleotide was called as belonging to an exon. 
This had the additional benefit of merging overlapping 
predicted exons . 

20 Finally, we applied a lower size threshold of 

75 contiguous nucleotides to each consensus exon. 

Each probe was completely sequenced on both 
strands prior to its use on a genome-derived single 
exon microarray; sequencing confirmed the exact 

25 chemical structure of each probe. An added benefit of 
sequencing is that it placed us in possession of a set 
of single base-incremented fragments of the sequenced 
nucleic acid, starting from the sequencing primer 3 1 
OH. (Since the single exon probes were first obtained 

30 by PCR amplification from genomic DNA, we were of 

course additionally in possession of an even larger set 
of single base incremented fragments of each of the 
single exon probes, each fragment corresponding to an 
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extension product from one of the two amplification 
primers . ) 

Hybridization analysis was conducted 
essentially as set forth in Examples 1 and 2, with one 
5 modification . 

In Examples 1 and 2, we used a pool of 10 
tissues/cell types as control. We have since observed 
that every probe that demonstrates expression in the 
control pool can readily be shown to be expressed in 

10 HeLa cells, and have used HeLa as the source of control 
message in the more recent experiments. 

In the analysis of hybridization results, the 
uniform absolute signal intensity threshold used in 
Examples 1 and 2 to identify signals large enough to be 

15 considered biologically significant (0.5, representing 
a level roughly 10 times greater than the average of 
all E. coli control spots on a first iteration chip) 
was replaced with a statistical threshold determined 
for each channel and each hybridization as follows. 

20 Starting typically with 32 E. coli sequences, 

spotted in duplicate (left and right side) for a total 
of 64 control spots per microarray, control spots were 
eliminated if we observed more than a five-fold 
difference between the left and right side raw 

25 (unnormalized) signals for the probe. 

The median of the normalized signal from the 
remaining control spots was calculated (see infra for 
normalization routine) . 

Control spots were eliminated as outliers if 

30 they had signal intensity greater than the median of 
the normalized signals plus 2.4 (where 2.4 is roughly 
12 times the observed standard deviation of control 
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spot populations) and normalization was performed as 
set forth below. 

The mean and standard deviation of the 
normalized signal intensity from the remaining control 
5 spots were calculated, and the mean plus three standard 
deviations of the controls was then applied as a 
minimum intensity threshold for the particular 
hybridization experiment, giving a 99% confidence that 
expression is significant. 

10 Signal normalization was accomplished as 

follows. For each hybridization (each microarray, 
separately for each of the two colors), the median 
value of all of the spots was determined. For each 
probe, the normalized signal value is the arithmetic 

15 mean of the probe ? s duplicate intensities (each DNA 

probe, including controls, is spotted twice per slide) 
divided by the population median. 

Using this threshold, we identified over 
15,000 single exon probes that produce significant 

20 signal in one or more of ten tested tissues/cell types. 
The exact structures of these single exon probes are 
clearly presented in the SEQUENCE LISTINGS included in 
commonly owned and copending U.S. provisional 
application nos . 60/207,456, filed May 26, 2000; 

25 60/234,687, filed September 21, 2000; 60/236,359, filed 
September 27, 2000; in commonly owned and copending 
U.K. patent application no. 24263.6, filed October 4, 
2000; and in commonly owned and copending PCT 
applications filed January 29, 2001 (attorney docket 

30 nos. PB 0004 WO 1, for "Human genome-derived single 
exon nucleic acid probes useful for analysis of gene 
expression in human heart"; PB 0004 WO 2, for "Human 
genome-derived single exon nucleic acid probes useful 
for analysis of gene expression in human brain"; PB 
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0004 WO 3, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in human adult liver"; PB 0004 WO 4, for 
"Human genome-derived single exon nucleic acid probes 
5 useful for analysis of gene expression in human fetal 
liver"; PB 0004 WO 5, for "Human genome-derived single 
exon nucleic acid probes useful for analysis of gene 
expression in human lung"; PB 0004 WO 6, "Human genome- 
derived single exon nucleic acid probes useful for 

10 analysis of gene expression in human bone marrow"; 
PB 0004 WO 7, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in human placenta"; PB 0004 WO 8, for "Human 
genome-derived single exon nucleic acid probes useful 

15 for analysis of gene expression in BT 474 cells"; 

PB 0004 WO 9, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in HBL 100 cells"; PB 0004 WO 10, for "Human 
genome-derived single exon nucleic acid probes useful 

20 for analysis of gene expression in Hela cells"), the 
disclosures of which are incorporated herein by 
reference in their entireties. 

We also predicted the sequence of the ORF 
within the exon of each of the probes, where ORF was 

25 defined as that portion of an exon that can be 
translated in its entirety into a sequence of 
contiguous amino acids. 

To predict the ORF, we first looked for 
consensus as between any two or more of the four gene 

30 prediction programs. Consensus was required in two 
parameters: (1) as with prediction of the exon, each 
nucleotide must have been identified by two or more 
programs as falling within an exon; and, additionally, 
(2) the programs relied upon to establish that 
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consensus must have agreed on the frame. Presence of a 
stop codon disqualified the predicted ORF. ORFs 
shorter than 50 nt were also disregarded. 

Absent consensus as to nucleotide and frame, 
5 each of the six frames of the predicted exon were 

examined individually for stop codons and the longest 
open reading frame of at least 51 nt selected as the 
exon's likely ORF. Certain of the exons have no ORF as 
defined by either set of criteria. 

10 We then translated the predicted ORFs using 

the standard genetic code. 

The exact structures of these single exon 
probes are clearly presented in the SEQUENCE LISTINGS 
included in commonly owned and copending U.S. 

15 provisional application nos. 60/207,456, filed May 26, 
2000; 60/234,687, filed September 21, 2000; 60/236,359, 
filed September 27, 2000; in commonly owned and 
copending U.K. patent application no. 24263.6, filed 
October 4, 2 000; and in commonly owned and copending 

20 PCT applications filed January 29, 2001 (attorney 

docket nos. PB 0004 WO 1, for "Human genome-derived 
single exon nucleic acid probes useful for analysis of 
gene expression in human heart"; PB 0004 WO 2, for 
"Human genome-derived single exon nucleic acid probes 

25 useful for analysis of gene expression in human brain"; 
PB 0004 WO 3, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in human adult liver"; PB 0004 WO 4, for 
"Human genome-derived single exon nucleic acid probes 

30 useful for analysis of gene expression in human fetal 
liver"; PB 0004 WO 5, for "Human genome-derived single 
exon nucleic acid probes useful for analysis of gene 
expression in human lung"; PB 0004 WO 6, "Human genome- 
derived single exon nucleic acid probes useful for 
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analysis of gene expression in human bone marrow"; 
PB 0004 WO 7, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in human placenta"; PB 0004 WO 8, for "Human 
5 genome-derived single exon nucleic acid probes useful 
for analysis of gene expression in BT 474 cells"; 
PB 0004 WO 9, for "Human genome-derived single exon 
nucleic acid probes useful for analysis of gene 
expression in HBL 100 cells"; PB 0004 WO 10, for "Human 

10 genome-derived single exon nucleic acid probes useful 
for analysis of gene expression in Hela cells")/ the 
disclosures of which are incorporated herein by 
reference in their entireties. 

The sequence of each of the probes, exons, 

15 and ORF-encoded peptides was used as a query to 

identify the most similar sequence in each of dbEST, 
GenBank NR, and SWISSPROT. The query programs used 
were BLAST (nucleic acid sequence query of dbEST and 
NR) , BLASTX (nucleic acid sequence query of SWISSPROT) , 

2 0 T BLASTX (peptide sequence query of dbEST and NR) , and 
BLASTP (peptide sequence query of SWISSPROT) . Because 
the query sequences are themselves derived from genomic 
sequence in GenBank, only nongenomic hits from NR were 
scored. 

25 The attached SEQUENCE LISTINGS in our 

commonly owned and copending applications report, for 
each SEQ ID NO:, the accession number of the entry from 
each of the three queried databases that gave the 
highest absolute expect ("E") value (the "top hit"), 

30 along with the "E" value itself. The SEQUENCE LISTING 
is incorporated herein by reference in its entirety. 



- 100 - 



All patents, patent publications, and other 
published references mentioned herein are hereby 
incorporated by reference in their entireties as if 
each had been individually and specifically 
5 incorporated by reference herein. While preferred 

illustrative embodiments of the present invention are 
described, it will be apparent to one skilled in the 
art that various changes and modifications may be made 
therein without departing from the invention, and it is 
10 intended in the appended claims to cover all such 

changes, modifications and equivalents that fall within 
the true spirit and scope of the invention. 



